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Abstract. The recent "breakdown criterion" result | 17| of S. Klainerman and 
I. Rodnianski stated roughly that an Einstein-vacuum spacetime, given as a 
CMC foliation, can be further extended in time if the second fundamental form 
and the derivative of the lapse of the foliation are uniformly bounded. This 
theorem and its proof were extended to Einstein-scalar and Einstein-Maxwell 
spacetimes in the thesis 1211 . In this paper, we state the main results of |21l , 
and we summarize and discuss their proofs. In particular, we will discuss 
the various issues resulting from nontrivial Ricci curvature and the coupling 
between the Einstein and the field equations. 



1. Introduction 

The general breakdown/continuation problem for PDE is the following: 

Under what conditions can an existing local solution of an evolution equa- 
tion on a finite interval [Tq,T) be further continued past T? 

This can be equivalently posed as the breakdown of such a solution at a finite time 
implying the violation of such conditions. The determination of such breakdown 
criteria can be a potentially useful step toward characterizing the blowup of solu- 
tions. Furthermore, in some instances, such breakdown conditions can be critical 
tools for proving global existence results. 

In this paper, we will consider this breakdown problem for the Einstein-scalar 
(E-S) and Einstein-Maxwell (E-M) equations in the CMC gauge. The results and 
proofs described in this paper are extensions of those in [T7] , which established the 
analogous breakdown criterion for the Einstein-vacuum (E-V) case. The matter 
held present in these nonvacuum cases presents additional issues to be addressed. 
The full details of this work can be found in the original report |21) . 

1.1. Classical Results. Consider, as a model example, the initial value problem 
for the following nonlinear wave equation on R 1+3 : 

(1) □<£ = (d t (f>) 2 , <p\ t=0 = 0o, <M t=0 = 0i. 

From classical theory, cf. [5D1 Thm. 22], we have the following results: 

• Given initial data (</> ,</>i) e Z s = H S (R 3 ) x H S -\R 3 ), where s > 5/2, 
then a unique solution to ([1]) exists in the space 

X. tT > = C ([0, T'] ; H s (R 3 )) n C 1 ([0, T'} ; H s - 1 (R 3 )) 

for sufficiently small T" > depending on the Z s -norm of the initial data. 

• The maximal time of existence T, i.e., the supremum of such r's for which 
a solution exists to time r, depends on the Z s -norm of the initial data. 
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A continuation result follows naturally as a companion to local well-posedness: 
Suppose a solution (j> on a finite time interval [0, T) to ([I]) is in X Sj t> for 
every < T' < T, where s > 5/2. If <fi also satisfies the criterion 

( 2 ) II^IIl°=([0,T)xR3) < °°, 

then 4> can be extended past time T as a solution of (HJ . Moreover, this 
extension is an element of X s ^+e for small e > 0. 
In terms of breakdown, this can be equivalently stated as follows: 

Suppose a solution </> on a finite time interval [0, T) to (|TJ) is in X s ,t' for 
every < T" < T, where s > 5/2. Then, if 4> breaks down at time T, i.e., 
if 4> cannot be extended as above, then d(f) $ L°°([0,T) x R 3 ). 
The main idea behind the result is the following observation: if ^ holds, then we 
can uniformly bound the £P(R 3 ) x £P _1 (R 3 )-norms of (<j>,dt<j>) on each timeslice 
{t} x R 3 , where < t < T. Therefore, we can apply the previous local well- 
posedness result to generate local solutions existing for a fixed time e > with each 
of the above cross-sections as the initial data. By uniqueness, we can patch these 
local solutions into a solution which exists on the interval [0, T + e). 

Remark. The condition @ is certainly not optimal, since slightly weaker iterated 
norm conditions are also known to be sufficient. 

A multitude of breakdown results have been established for other evolution equa- 
tions. For example, consider the incompressible 3-dimensional Euler equations 

(3) d t u + u ■ Vu + Vp = 0, V • u = 0, 

where u : R 1+3 — ^ R 3 represents the velocity and p : R 1+3 — > M represents the 
pressure. In addition, define the vorticity of u to be the curl u> = V x u of u. 

A well-known result of Beale, Kato, and Majda in [3] established that if a local 
finite-time solution has its vorticity bounded in the i(L^°-norm, then the solution 
can be extended further in time. An important point here is that unlike the 
condition J2]), we need not bound all components of the derivative of the solution. 
The proof is in principle like that of the nonlinear wave equation; we use this L\U^- 
bound on the vorticity in order to derive uniform energy bounds related to the local 
well-posedness theory of the equations ([3|) . 

Another example of a breakdown condition lies in the paper [S] of Eardley and 
Moncrief on the Yang-Mills equation in R 1+3 . In this setting, a sufficient continu- 
ation criterion is an L°°-bound on the Yang-Mills curvature F. Using the standard 
representation formula for the wave equation, however, one can demonstrate that 
such a uniform bound always holds. The immediate consequence, then, is a global 
existence result. Furthermore, Chrusciel and Shatah, in [5], generalized this result 
to globally hyperbolic (1 + 3)-dimensional Lorentzian manifolds using mostly the 
same principles, but applying instead the representation formula of [llj . 

1.2. Results in General Relativity. In general relativity, a number of break- 
down results have been established for the E-V equations. For example, by solving 
the equations in the standard fashion by imposing favorable gauge conditions, such 
as wave coordinates, one can show that given a solution in which an L°°-bound 

^The L\ LJ°-norm of u is the Z^-norm of the function 1 1— ^ |]w(t)|| £ oo m3) ■ 
2 For simplicity, we neglect the Higgs field, which was also discussed in [9]. 



BREAKDOWN CRITERIA 



3 



holds for d<3, where © is the spacetime metric and d refers to the gauge coordinate 
derivatives, then the solution can be further continued. Examples of such results 
include [6] and, more recently, [2J. Although such a condition is quite analogous 
to the model case of the nonlinear wave equation, it is also non-geometric, as it 
depends on the specific choice of coordinates. Moreover, the condition requires 
bounds on all components of the derivative of the metric. 

Another more geometric breakdown result for the E-V equations was given by M. 
Anderson in [T] . Here, the continuation criterion is an L°°-bound on the curvature of 
the spacetime. While this is clearly geometric, it does have the added disadvantage 
of depending on essentially two derivatives of the spacetime metric. 

The next point of discussion is the improved breakdown result of S. Klainerman 
and I. Rodnianski for E-V spacetimes, presented in [17]. Their main result, stated 
in [17[ Thm. 1.1], can be summarized as follows: 

Theorem 1. Suppose (M,g) is an E-V spacetime, given as a CMC foliation 
M= (J E T , t <h< 0, 

t <T<t t 

where each S r is a compact spacelike hyper surface of M satisfying the constant 
mean curvature condition tr k = r, and where k denotes the second fundamental 
form of E r in M . In addition, let n denote the lapse of the S T 's, and assume the 
following breakdown criterion holds: 

(4) \\ k L°°[M) + l|V(logra)|| ioo(M) < oo. 

Then, the spacetime (M,g) can be extended as an E-V spacetime in the CMC gauge 
to some time ti + e for some e > 0. 

We make note of the following important features of Theorem [T] 

• This result is more geometric than the previous coordinate breakdown con- 
dition, whose statement required the choice of an entire coordinate system. 
TheoremlU on the other hand, relies only on the (constant mean curvature) 
time foliation to state the breakdown criterion. 

• Both k and V(logn) in Q reside at the level of one derivative of g. Also, 
these quantities do not represent all components of one derivative of g. 

The proof of Theorem [1] is a complicated affair which in its totality spans sev- 
eral papers: [HI [13l [TH [15j [16j [17l HH [25]. One important point of the proof is 
that although the associated local well-posedness problem applies to objects on the 
timeslices, much of the technical work revolves around spacetime objects, in par- 
ticular the spacetime curvature. The relevant spacetime and timeslice objects can 
then be related using standard elliptic estimates. 

Using mostly the same set of ideas as in [T7], D. Parlongue, in [19] , proved 
an analogous breakdown criterion in the setting of a maximal time foliation with 
asymptotically flat timeslices. In addition, the breakdown condition was weakened 
from an L°°-bound to an iterated L^L^-bound: 

\\ k hjL^(M) + ll V ( lo g n )IL?L~(M) < °°' 

Furthermore, there have been additional results by Q. Wang toward an improved 
£t L^-criterion in the original CMC setting; see [25 1 127 ] , 
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1.3. Nonvacuum Spacetimes. The main question posed here is the following: 

Do breakdown results analogous to that of Theorem [7] hold for Einstein- 
scalar and Einstein- Maxwell spacetimes? 

In this paper, we will answer this question affirmatively. 

The breakdown criteria in these nonvacuum cases remain largely the same as in 
the vacuum case. We retain the uniform bounds on the time foliation quantities k 
and V(logrt). However, we must also impose a uniform bound on the now nontrivial 
matter field: the scalar field in the E-S setting, or the spacetime Maxwell 2-form 
in the E-M setting. This will be the only additional condition. 

In the vacuum case of |17j . the breakdown criterion in conjunction with the CMC 
gauge imply a wide range of priori controls with respect to the time foliation, all 
of which are essential for proving Theorem [TJ These include, for example, uniform 
L 2 -bounds for curvature along the timeslices, as well as uniform bounds for Sobolev 
constants on the timeslices. In the E-S and E-M cases, we will also need uniform 
bounds on the matter fields in order to achieve the a priori controls mentioned 
above. This will justify the modified breakdown criterion which we shall adopt. 

We note that many of the elements which made the vacuum case so difficult will 
also make our nonvacuum cases similarly demanding. For instance, both |17j and 
the main result here are essentially "large data" results. Therefore, the only a priori 
"energy estimates" we will have are lower-order L 2 -bounds on the curvature and 
the matter field. This greatly complicates the process behind controlling the local 
null geometry, since both the null injectivity radius and the Ricci coefficients will 
need to be controlled by these L 2 -quantities. In particular, this will necessitate the 
use of the geometric Littlewood-Paley theory and the associated Besov estimates 
developed in [13] and further discussed in [21] Sec. 2.2]. 

The addition of matter fields also introduces a number of new difficulties: 

• The spacetime Ricci curvature is now nontrivial. 

• We must also deal with the coupling between the curvature and the matter 
field. An unavoidable consequence of this is that both the curvature and 
the matter field must be estimated concurrently. 

• In the E-M case, there exist first-order terms in the wave equations satisfied 
by the curvature and Maxwell field which were not present in [17] . These 
terms cannot be fully treated using only the techniques of [17] and its sub- 
sidiary papers. In particular, we will require a generalized representation 
formula for tensor wave equations; see [21] Ch. 5] and [22]. 

Another unfortunate reality is that the assumption of an E-V spacetime was 
pervasive throughout the entire proof of Theorem [T] As a consequence, most of 
the elements of the proof in the E-V case must in principle be redone. In particular, 
this includes controlling the local null geometry. In both the vacuum case of [17] 
and the nonvacuum cases presented here, the work on the null geometry comprises 
the most lengthy and technically involved portion of the overall proof. We will 
omit this portion of the proof from this paper, since it would more than double the 
length of this text. For details in this area, see [21], as well as its predecessors in 
the vacuum setting: Q2] Q3] QU QU [24] [25] . 

The following statement summarizes the main theorem of this paper: 



■^Notable exceptions include the geometric Littlewood-Paley theory of I13| and the Kirchhoff- 
Sobolev parametrix of [15] (and its generalization in [211 Ch. 5] and |22|). 
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Theorem 2. Suppose (M,g, $) is an E-S or E-M spacetime, given by 
M = (J S T) t < t x < 0, 

t a <T<t! 

where each S r is a compact spacelike hypersurface of M satisfying the constant 
mean curvature condition tr k = r , and where $ denotes the matter field: 

• In the E-S case, $ denotes the scalar field 4>. 

• In the E-M case, denotes the Maxwell 2- form F . 

• Let k and n be as before in the statement of Theorem [7J 
In addition, let the quantity $ denote the following: 

• The spacetime covariant differential D(f> of <f>, in the E-S case. 

• The Maxwell 2-form F, in the E-M case. 
Assume the following breakdown criterion holds: 

(5) ll*lli-(M) + l|V(Iogn)|| ioo(M) + ||ff|| i » (M) <oo. 

Then, (M,g,<&) can be extended as an E-S or E-M (resp.) spacetime in the CMC 
gauge to some time tx + e for some e > 0. 

A more technically precise version is stated as Theorem [51 after a sufficient 
amount of background and notations have been developed. 

1.4. Extensions and Open Problems. Finally, we discuss some possible exten- 
sions of Theorem (5J as well as some related open problems. The most immediate 
question is whether Theorem [5] can be similarly adapted to other related classes 
of nonvacuum spacetimes, such as Einstein-Klein- Gordon and Einstein- Yang-Mills 
spacetimes. The answer is affirmative for the above two spacetimes in this same 
CMC gauge framework. The class of techniques presented in this text can be ap- 
plied directly to the Einstein-Klein-Gordon case. This is also mostly true for the 
Einstein- Yang-Mills setting; however, in this case, we will also require vector bundle 
generalizations of many of the tools used here to handle the Yang-Mills curvature 
both in a covariant and in a gauge-invariant fashion. 

Whether analogous breakdown criteria can be stated and proved for less simi- 
lar classes of nonvacuum spacetimes, such as the Einstein-Euler and the Einstein- 
Vlasov models, remains an open question. For such matter fields, one would likely 
require other tools besides the representation formula for tensor wave equations in 
order to derive higher-order energy inequalities. 

Another variation of Theorem[2]which should be possible is the case of a maximal 
foliation with asymptotically flat timeslices for various nonvacuum spacetimes. The 
adaptations to the theorem statement and proof should then be analogous to [19] . 
Note that one must also amend the energy norms for the curvature and the matter 
fields in this setting along the lines of [19] . 

One more potential direction of investigation involves the possible weakening of 
the L°°-bounds in the breakdown criterion to slightly weaker iterated bounds, such 
as the L^L^ bounds found in [T5]. A further improvement to Lji^-bounds along 
the lines of [551 HZ] is also expected to be achievable. 

A more difficult open question related to the topic at hand is the well-known 
CMC conjecture. In this case, the desired conclusion is that the spacetimes of 

4 In particular, this includes the vector bundle extension of the generalized Kirchhoff-Sobolev 
parametrix; see |2 1 1 Sec. 5.3] and [221 Sec. 5]. 



6 



A RICK SHAO 



Theorems Q] and [2] can in fact be continued up to time 0. Finally, another related 
conjecture is the bounded L 2 -curvature conjecture, which in principle states that 
a local well-posed theory exists based on control of the L 2 -norms of the curvature. 
There has been significant recent progress in resolving this question. 

Acknowledgements. The author wishes to thank Professor Sergiu Klainerman 
for suggesting this problem and for hours of discussions throughout the preparation 
of this work. Thanks also extends to Professor Igor Rodnianski for his insights. In 
addition, the author thanks Qian Wang for helpful technical discussions. 

2. Tensorial Notations 

In this section, we construct some notations for the various types of tensor fields 
we will encounter. The main objective of these notations is to highlight the co- 
variant structures present in our setting. @ Although such structures were also 
pertinent to the analogous vacuum problem of [T7], as well as other related works 
(for example, [3 [T^l HHDi they were described more implicitly in those texts. We 
develop the new and more explicit notations here in the hopes of promoting a more 
conceptual outlook on the objects of our analysis. The conventions used here will 
be an abridged version of that of [21] . 

2.1. Tensor Fields. Let M denote an arbitrary smooth manifold. First, we set 
the notations for standard tensor bundles over M: 

• Let TJM denote the rank-(r, s) tensor bundle of M. 

• Let S S M denote the bundle of all fully symmetric elements of T®M. 

• Let A S M denote the bundle of all fully antisymmetric elements of T°M, 
i.e., the exterior bundle of degree s on M. 

In general, for a smooth vector bundle V over M , wc dchnc the following: 

• Given p £ M, let V p denote the fiber of V at p. The lone exception is for 
tangent spaces, for which we use the more standard notation 



• Let TV denote the smooth sections of V. For example, TTJM and TA S arc 
the spaces of rank-(r, s) tensor fields and s-forms on M, respectively. 
For convenience, we also define the abbreviations 



i.e., the space of vector fields on M, the space of 1-forms on M, and the space of 
all tensor fields on M, respectively. 

When convenient, we will adopt standard index notation to represent tensor 
fields with respect to local frames and dual coframes. In many instances, these 
frames satisfy additional conditions, e.g., coordinate, orthonormal, or null frames. 
In accordance with Einstein summation notation, indices repeated in both super- 
script and subscript represent a summation of components, i.e., a contraction. In 
addition, we will use capital letters to denote collections of indices. Moreover, 
repeated capital letters signify summations over all represented indices. 



'By covariant structures, we mean vector bundles with metrics and compatible connections. 




(6) 



X(M) = TT$M, X* (M) = rT 1 °M, TTM = |J TT r s M 



r,s>0 
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Furthermore, when writing in index notation, we let 2t Q| g [•] denote an unnormal- 
ized antisymmetrization of the indices a and j3 of the indexed quantity within the 
bracket. For example, if S £ TT^M and T £ TT°M, then 

(7) 2l^ 7 [S a pT 7 s] = S a pT~[S - S a ~ l Tps — 2S a [pT^ 5 . 

Although the standard bracket notation on the right-hand side of ([7]) achieves 
the same effect, the 21- notation will be useful for larger expressions with multiple 
antisymmetries, in which the bracket notation would be confusing or ambiguous. 

In the case that M is Riemannian or Lorentzian, with metric g £ TS 2 M, then 
contravariant and covariant components are equivalent, and we will often write 
T r+S M for TJM. Moreover, given a field T £ TT r M, we define R[T] £ TT r+2 M 
to be the covariant derivative commutator 

(8) R a p [T]j = D a pT[ - Dp a Ti. 

We will also follow these standard traditions within index notation: 

• We denote both the Ricci and scalar curvatures of (M,g) by R, but with 
the appropriate number of indices. 

• If M is oriented, then the volume form is denoted e. 

Finally, if (M, g) is Riemannian, we can define natural tensor norms by 

(9) |-| : YTM -» C°° (M) , |*| = (V*,)^ . 

We can then define the obvious L p -norms, I < p < oo, for such tensor fields. 

2.2. Foliations. Next, we briefly discuss a rather general situation of a I-parameter 
foliation of Riemannian manifolds. Let (M, g) denote a Lorentzian manifold, and 
let N denote a smooth submanifold of M. Let / £ C°°(N), mapping onto a possibly 
infinite open interval X, with df nonvanishing. For each x £ I, we assume 

S, = {?e M \f( P ) = x} 

is a Riemannian submanifold of M . In other words, we can write 

N= |J 

as a 1-parameter foliation of Riemannian submanifolds. 

In future sections, this situation will arise in the following instances: 

• Time foliations of spacetime: In this case, N = M, and / is a time function 
on M which assigns to each a point a time value; see Sec. 13.11 

• Spherical foliations of regular null cones: In this case, N is a smooth portion 
of a past null cone in M, and / foliates N into spherical cross-sections; see 
Sec. 14.21 In particular, N is null and hence is not pseudo-Riemannian. 

Since it will be convenient to have some common notations for both cases, we 
commit to this abstract development here. 

We begin by defining the following bundles over N: 

• A tensor w £ (T r M) q at q £ Y> x is horizontal iff w is tangent to E x . 

• We denote by T_ k N the horizontal bundle over N of all horizontal tensors 
of total rank k at every q £ N. 

• We denote by T l N the extrinsic bundle over N of all tensors in M of total 
rank I at every q £ N, i.e., the restriction of T l M to N. 



'For example, the Ricmann curvature _R of (M, g) is an element of FT 4 M. 
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• The mixed bundle T l T k N over N is defined to be the tensor product bundle 

(10) T k T l N — T k N ® T l N. 
For convenience, we adopt the notations 

(11) TTN = \J rr fe iV, TTN = y TT l N, TTTN = [J TT k T l N : 

k>0 l>0 k,l>0 

i.e., the spaces of horizontal, extrinsic, and mixed tensor fields on N. We also define 

(12) X (AO = TT}N, X (N) = TT X N. 

Moreover, we adopt the following general indexing conventions: 

• Horizontal indices will be denoted using Latin letters. 

• Extrinsic indices will be denoted using Greek letters. 

• Collections of extrinsic indices will be denoted using capital Latin letters. 
For example, the pullback to N of g is an element of TT 2 N, while the induced 

metrics 7 on the E^'s is in rT 2 iV. Note that 7 induces positive-definite bundle 
metrics (•, •) on the horizontal bundles T fc 7V via full metric contraction. Similarly, 
the restriction of g to N induces bundle metrics (•, •) on the extrinsic bundles T l N. 
These in turn naturally induce "product" bundle metrics on T_ k T l N. 

Next, let D denote the Levi-Civita connection on (M,g). Recall from standard 
theory that D induces a connection D on the extrinsic bundles T l N. Moreover, 
Dxg = by definition, hence D is compatible with the extrinsic bundle metrics, [j 

We can also naturally define a connection V on the horizontal bundles via pro- 
jections. In particular, given X £ X(N), then V satisfies the following properties: 

• If / e C°°(N), then V x / = Xf, as usual. 

• If Y £ X(N), then VxY is the projection onto the E^'s of DxY. 

• If A £ TT k N is fully covariant, and Y x , . . . , Y k £ X(N), then 

(13) V x A(Y 1 ,...,Y k ) = X[A(Y 1 ,...,Y k )}-A(VxY 1 ,Y 2 ,...,Y k ) 

- A{Y 1 ,...,Y k _ 1 ,V x Y k ). 

This definition of V generalizes the usual Levi-Civita connections on the E^'s to also 
include non- horizontal derivatives tangent to N. Since Vx7 = for any X £ X(N), 
then V remains compatible with the horizontal bundle metrics. 

We can also define mixed connections V on the mixed bundles using the above 
connections V and D. On TJ^_ k T l N, we define V to be the unique connection 
satisfying the following Leibniz identity for decomposable fields: 

Vx (A® B) = V X A ®B + A® D X B, X £ X (N) , A £ TT k N, B £ TT l N. 

In other words, V behaves like V and D on horizontal and extrinsic components, 
respectively. Clearly, V is compatible with the above mixed bundle metrics. 

Remark. More explicitly, if A £ TT k T l N; X £ X(N); Y u ...,Y k £ X(N); and 
Z\, . . . , Zi £ X(N); then V xA(Yi, . . . , Y k ; Z\, . . . , Zi) is assigned the value 

(14) X [A {Yi, .. ..)),:/, Z t )] - A (V X Y U Y 2 , ....)),: Z. ... Z, : ... 

- A (y x , . . . , r fc _i, V x Y k ; Z X ,...,Z X )-A (Y u . . . , Y k ;D x Z 1 ,Z 2 , ...Z x ) 
-...-A(Y 1 ,...,Y k ;Z 1 ,...,Z l _ 1 ,D x Z l ). 



^Technically, by Dxg, we mean Dj acting on the restriction of g to N. 
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We also set the following notations: 

• For A e TT k N, define VA <E TT k+1 N to map X € X(N) to V X A. 

• For A £ TT k T l N, define g TT k+1 T l N to map X € X{N) to VxA 

• Horizontal and mixed Laplacians are defined in the usual fashion: 

(15) A = 7 ab V ab , A = 7 ab V a6 . 

For further details involving the above constructions, see [2TJ Sec. 1.2]. 

Remark. In contrast, the wave operator g a ^D a p on (M,g) is denoted □. 

The fundamental properties are the Leibniz rules satisfied by the connections and 
the compatibility between the connections and metrics. These justify integration 
by parts operations involving V-derivatives, seen in [TSJ (T7] [3TJ [35] and later in this 
paper. This was used implicitly in [H1[T7] and was discussed in detail in [2T1 [33] . 

2.3. Normal Transport. Assume the same foliation setting as before. For each 
point p €E Tt x C N, there is a unique direction in N which is normal to S x . Thus, 
we can define the "normalized normal" vector field Z £ X(N) such that Z points 
in the direction in N normal to the S x 's and satisfies Zf = 1. 

We can naturally define quantities via transport along the integral curves of 
Z. For instance, given a locally defined function on some E y , we can define a 
corresponding local function on another E x via the diffeomorphisms induced by 
the above transport. In particular, given a coordinate system (y 1 ,?/ 2 ) in E a , we 
can define a transported coordinate system on the other E^'s. 

Similarly, we can transport horizontal tensors along the integral curves of Z. 
This induces the standard definition of "Lie derivatives" of horizontal tensor fields. 
More explicitly, for any A e TT_N, its "normal Lie derivative" CfA is defined 



(16) C f A\ p = Jim v f -, pe^CJV, 

where f x +s,x is the flow from Ti x+ $ to T, x via the integral curves of Z , and where 
fx+Sx denotes the "push- forward" of tensors through f x +s,x- 

Remark. The development described here differs from that of [171 121] . which ex- 
pressed similar notions entirely in terms of transported coordinate systems. Notice 
that the normal Lie derivatives of (|16[) . which are tensorial and invariant, coincide 
with the f -coordinate derivatives for such transported coordinate systems. 

Finally, we define the analogous "normalized normal" covariant derivatives: for 
fields A € TTN and B G TTTN, we define 

(17) X7 S A = V Z A, V f B = V z B. 

2.4. The Einstein Equations. Assume now that (M,g) is a (1 + 3)-dimensional 
connected, oriented, and time-oriented Lorentzian manifold, with Levi-Civita con- 
nection D and Riemann curvature R. In addition, let f denote a collection of pre- 
scribed "matter fields" on M, satisfying the equations of their respective theories. 
The Einstein equations on M are given in index notation by 

(18) G a p = R a p — — ■ R ■ g a (j = Q a p, 
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where Q G TS 2 M denotes the energy-momentum tensor of the matter fields. Recall 
that the Ricci curvature can be expressed in terms of the matter fields: 

(19) R a p = Qap - ^gapg^Q^- 

Moreover, recall that Q is both symmetric and divergence-free. 

In this text, we will consider the following settings for matter fields: 

• Vacuum: There are no matter fields, so both Q and Ric vanish. 

• Scalar field: In this model case, f is given as a scalar field <fi G C°°[M), 
which satisfies the linear homogeneous covariant wave equation 

(20) Ug<j) = 0. 

The energy-momentum tensor and Ricci curvature are given by 

(21) Q aP = D a 4>Dp4> - -gapIPtDpfa R a0 = D a 4>D p 4>. 

• Maxwell field: In the model electromagnetic case, f is given as a Maxwell 
field F G TA 2 M, which satisfies the linear homogeneous Maxwell equations 

(22) D^F a p = 0, D a F M + D p F ia + D 7 F aP = 0. 

The energy-momentum tensor and Ricci curvature are given by 

(23) Q a p — R a p — Fa^F 13 ^ - -gapF^F^is- 

We remark that in the above settings, Q satisfies both the positive energy con- 
dition and the strong energy condition, i.e., Q(X, Y) > and Ric(X, Y) > for all 
future causal X,Y € X(M). Moreover, by direct calculations, we can compute 

(24) D a R a p 7S = 0, 

(25) D a R aPjS = -% s [D^D sp <f>] , 

(26) D a R aMS = -Fp^DpFys - %s F^D s F 0tl + ^g^F^DsF^ 
in the E-V, E-S, E-M settings, respectively. 

3. The CMC Breakdown and Cauchy Problems 

For the breakdown problem at hand, we are given an existing solution of the 
Einstein equations, i.e., an E-S or E-M spacetime "existing on a finite time interval" . 
Our goal will be to "extend the solution further in time" given that certain criteria 
hold on this existing solution. In this section, we will make precise the above 
informal expressions. This essentially requires smoothly assigning to each point of 
the spacetime a "time value" and then foliating the spacetime into "timeslices" , 
that is, hypersurfaces of constant time value. 

This objective will be achieved by defining global time functions and the various 
fundamental objects associated with its resulting foliation. In addition, for the 
breakdown problem, we will adopt the CMC gauge condition as well. 

Throughout the remainder of this paper, we will always let (M, g) denote a 
(I + 3)-dimensional connected, oriented, and time-oriented Lorentzian manifold, 
with Levi-Civita connection D and curvature R. 



8 See Q. 
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3.1. Time Foliations. We define a time function on M to be a map t G C°°(M) 
satisfying the following conditions: 

• The spacetime gradient of t is everywhere past timelike. 

• Every nonempty level set of t is a Cauchy hypersurface of M, i.e., every 
inextendible causal curve in M intersects each level set of t exactly once. 

From the first property, we see t is strictly increasing along all future directions. 
Moreover, the second property implies that (M, g) is globally hyperbolic. 
For r G R and an interval ICR, define 

S r = {z G M | t (z) = t} , S/ = {z G M | t (z) G /} . 

The family {S T } of spacelikc Cauchy hypersurfaces defines a time foliation of M, 

(27) M=|JX T) 

rex 

where 1 is an interval in R, and where t maps onto I. We will also impose the 
following assumptions: the interval X is finite, and the E r 's are compact. 

Remark. Any two slices £ Tl , E T2 , where tj. , tj G I, are m /aci diffeomorphic, 
since any p\ G S Tl can be canonically identified with the point pi G St 2 on the 
integral curve of grad t through p\ . 

We adopt the notations developed in the Section l2~2"l to describe the foliation (|2"7| , 
with N = M and f = t. In addition, we adopt the following indexing conventions: 

• Greek letters refer to all components in M, ranging from to 3. 

• Latin letters refer to horizontal components, ranging from 1 to 3. 

• For implicit index summations, repeating Greek indices are summed from 
to 3, while repeating Latin indices are summed from 1 to 3. 

Also, let 7, V, and 1Z denote the horizontal metrics, connections, and curvatures. 
Next, we define the following basic objects of interest: 

• Let n G C°° (M) denote the lapse function, given by 

(28) n= \g (grad t, grad t)\~* > 0. 

• Let T G X(M) denote the future unit normal to the E T 's, i.e., 

(29) T=-n-gradt. 

• Let k G TT 2 M be the future second fundamental form of the S T 's: 

(30) k{X,Y) = -g{D x T,Y), X,YeX(M). 

Recall that k is symmetric and can be decomposed into trace and traceless parts: 

1 
3 

Recall also that tr k corresponds to the mean curvature of the E T 's in M. 

We can relate spatial derivatives of 1Z and k to spacetime derivatives of R. First, 
the Gauss equation, expressed in index notation in terms of TZ and k, becomes 

(32) Rijlm — T^ijlm ^jl^im ~l~ ^il^jm- 

Similarly, the Codazzi equations and a direct calculation imply 

(33) V'fcy = R ai T a + V, (tr k) , V {l k 0]l = R a ujT a . 



(31) tr k = j ij kij G C°° (M) , k = k - - (tr jfe) 7 G TT 2 M. 
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In this setting, the "normalized normal" vector field Z is given by Z = nT. In 
addition, the quantities 7 and k satisfy the evolution equations 

(34) Ca = -2nk, 

(35) Ctktj = -Van + nR la][j T a T p - nk u k/, 
where C t denotes the normal Lie derivative of ([TBI) . 




In other words, we define h to act like g on the horizontal components, and we 
simply invert the sign of the normal timelike component. It is easy to see that h is 
Riemannian, and that h agrees with 7 for horizontal vector fields. From now on, 
tensor norms on M will be defined with respect to h. 

The following proposition states that the Cauchy-Schwarz inequality continues 
to hold with respect to contractions by g: 

Proposition 3. Let e TTM, and let $ • * g YTM denote a tensor field 

obtained by taking zero or more contractions and g- contractions of <E> <g> \P . Then, 
|$ . < |$||^r| ; where the constant depends on the number of contractions. 

Proof. See (2TJ Prop. 4.1]. □ 

Let 7r = ( T '-7r 6 TS 2 M denote the deformation tensor of T: 



(37) ■K{X,Y)=C T g{X,Y)=g{D x T,Y)+g{D Y T,X), X,YeX(M). 



By direct computation, we can relate tt to n and k: for any X, Y G X(M ), 



(38) tt (T, T) = 0, tt{X,Y) = -2k(X,Y), it (T, X) = X (log n) . 



In other words, 7r can be thought of as a spacetime object which contains precisely 
the same information as the horizontal objects k and V(logn). 

One reason for the importance of the deformation tensor is that it bounds a 
number of essential quantities related to the time foliation, the most important of 
which are demonstrated by (|38"1) . We also have from [2TJ Prop. 4.2] the bounds 



For instance, the following calculus estimate is a direct consequence of (|39)): 
Proposition 4. The following calculus estimate holds: 



Proof. Recalling the Leibniz and metric compatibility properties of the mixed co- 
variant differential V, along with Proposition [3l we compute 



(39) 



DT\ < \n\ , 



\Dh\ < [tt 



P > i, * g rrrs T . 




*)] 



*| p ~ 2 \Dh\ |*| 2 . 



The proof is completed by applying (|39[) . 



□ 



Lastly, we note the following coarea formula: 



'See Section |2. 31 
By A < B, we mean A < CB for some positive constant C. 
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Proposition 5. IfQ is an open subset of M and J is a subinterval of I, then 

I 4 = I (I n-Adr 
for any integrable <j) £ C°°(VL). 

Proof. See (2U Prop. 4.4]. □ 

3.2. The Main Theorem. In terms of describing the Einstein equations as a 
system of partial differential equations, the above prescription of a time function 
and foliation leaves us with still another degree of freedom: a gauge condition 
satisfied by the foliation. In this text, we will adopt the constant mean curvature, 
or CMC, gauge, as was done in [17] . 

To be much more precise, we say that the spacetime (M,g), along with a time 
function t e C°°(M) and negative real numbers to < ti < 0, satisfy the condition 
(CMC)t l i 0i t 1 iff the following conditions hold: 

• The total timespan I = t(M) is precisely the interval (to,ti). 

• The timeslices S r are compact. 

• The CMC condition holds, that is, trk — t everywhere on M. 
We are now ready to state the precise main theorem of this text: 

Theorem 6. Suppose [M, g, 4>) is an Einstein-scalar or an Einstein-Maxwell space- 
time, where $ denotes the matter field <f> £ C°°(M) or F £ TA 2 M, corresponding 
to the E-S and E-M cases, respectively. In addition, suppose the following hold: 

• The condition fCMCj f i0iil holds for (M , g) . 

• The following "breakdown criterion" holds for some constant Co > 0, 
(40) \\k\\ L ~ { M) + l]V(logn)]] ioo(M) + \\3\\ L ~ (M) < Co, 

where $ denotes either D(j) or F in the E-S and E-M cases, respectively. 
Then, (M, g,&) can be extended past time t\ as a CMC foliation. In other words, 
there is an Einstein- scalar or Einstein- Maxwell (resp.) spacetime (M*, g*, $*) sat- 
isfying CCMCji^j^t^t, where t* is a time function on M+, e > 0, and t\ + e < 0, 
such that the following statements hold: 

• There exists an isometric imbedding i from (M, g) into (M*,<jr*). 

• The maps t and t+ correspond with respect to i, i.e., for every to < r < t\, 

i (S r ) = {qe M* | ti, (q) =t}, i (M) = {q e M* \t <U (q) < h} . 

In particular, we have t = t* o t. 

• The matter fields $ and also correspond with respect to i, i.e., $ = i*$*. 

The goal of the remainder of the paper will be to prove Theorem |6] 

3.3. The Cauchy Problem. We now give an explicit formulation of the Cauchy 
problems for the E-S and E-M equations in the CMC gauge. We will then state a 
basic local well-posedness result for this problem. 

A 4-tuple (So, 7o, fcoj fo) will be called an "admissible initial data set" for the 
Einstein-scalar or Einstein-Maxwell equations iff the following hold: 

• (So,7o) is a 3-dimensional compact oriented Riemannian manifold. 

• The field k is an element of r5 2 S . 



fco is to be interpreted as a "second fundamental form" for Eq. 
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• The "mean curvature" tr ko has a constant value tq < on all of So. 

• The field fo corresponds to initial data for the given matter field. In the 
E-S case, this is expressed as a pair <po,(/)i <E C°°(So), while in the E-M 
case, this is given as a pair Eq, Hq £ X*(So). 

• Letting V and TZq denote the Levi-Civita connection and the curvature for 
So, respectively, then in the E-S case, the following constraints hold, 

(41) - ^Vito = V J ' (fc )y , |V</> o | 2 + 0i =fto- N 2 + r 2 , 
while in the E-M case, the following constraints hold, 

(42) - e/ (E ) J (H )t = V j (fco)y , \E Q \ 2 + \H Q \ 2 =K - \k \ 2 + r 2 . 

• In the E-M case, the "matter field" f solves the additional constraints 

(43) V 1 (E ) i ee 0, V 1 (Hq), ee 0. 

Remark. We make the following remarks about the above definition: 

• The symbols | • | in [|41jl. (|42| denote the jo-tensor norm. The volume form 
tijk in (H2J is that of (So, 70), with respect to a chosen orientation 0/S0. 

• Note that (|41 [) and (|42p correspond to the E-S and E-M (resp.) constraint 
equations in the CMC gauge. Similarly, (|43[) corresponds to the standard 
constraints for the Maxwell equations. 

• If <fi and F are the desired "spacetime fields" for which we wish to solve, 
then 4>q and 4>\ correspond to the values <p and Ct4> on So, while Eq and 
Hq correspond to the electromagnetic decomposition of F on So . 

Now, if we are given an admissible initial data set (So, 70, ^Oj fo) m either the 
E-S or the E-M setting, then the goal will be to solve for a triple (M, g, <!>), along 
with maps i : So — » M and t g C°°(M), where the following hold: 

• (M,g) is a (1 + 3)-dimensional globally hyperbolic Lorentzian manifold. 

• There exists a time function t on M, along with constants tQ, t\ £ M, where 

< T o < ti < 0, such that (M,g) satisfies the condition (CMC) t t0itl . 

• The map i is an isometric imbedding of So into M, and i(So) is precisely 
the level set S T[) of t. 

• The element $ represents the matter field on M: in the E-S setting, then 
* = (j) € C°°(M), while in the E-M setting, we have $ = F G TA 2 M. 

• Both the Einstein equations (TT51) and the appropriate field equations for $ 
(either ^ or are satisfied on (M,g). 

• The field kg corresponds to the future second fundamental form of S To . In 
other words, fco coincides with the pullback i*k. 

• The field fo corresponds with the restriction of $ to So- To be more precise, 
in the E-S setting, this means 

4> Q = i*4>, 0i = \* c t 4>. 

In the E-M setting, if E, H is the electromagnetic decomposition of F, then 

Eq = i* E, Hq — i H. 

If all the above conditions hold, then (M,g, $), along with i and t, will be called a 
solution of the (Einstein-scalar or Einstein-Maxwell) CMC Cauchy problem corre- 
sponding to the initial data set (Sq, 70, ^0, fo)- 
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Our next task is to state a local well-posedness theorem for both of the above 
CMC Cauchy problems. This will be a straightforward modification of the vacuum 
analogue stated in [T71 Prop. 6.1]. 

Theorem 7. Let (So, 70, fccb fo) be an admissible initial data set for the E-S or 
E-M equations, and let To = tr fco < 0. Then, there exists a solution (M, g, 
of the E-S /E-M (resp.) CMC Cauchy problem corresponding to the above initial 
data set, with (M,g) satisfying (CtAC)t t a ,t\ f or some time function t on M and 
to < To < t\ < 0. Furthermore, the solution is unique up to isometric imbedding, 
and the time of existence t\ — tq depends continuously on the following parameters: 

• The initial mean curvature tq. 

• The diameter and injectivity radius of So ■ 

• The following Sobolev norm for fco : 

£0 = l|fco|| L 4 (Eo) + l|Vfc || L 2 (Eo) + ||V 2 fc || L2(So) + ||V 3 fc || i2(So) • 

• The following Sobolev norm for the curvature TZq of So •' 

^0 = ll^o|| L 2(So) + ll V7 ^o|| L 2 (So) + ||V 2 ft || i2(So) . 

• In the E-S case, the following Sobolev norm for fo = (</>o, 4>i) : 

So = l|V^o|| L 4 (So) + ||V 2 ^ || i2(So) + ||V 3 ^ || i2(Eo) + ||V 4 0o|| L2(s) 

+ Hi\\l H z o) + HV0i|| i2(So) + ||v 2 0i|| L2(So) + ||v 3 ^|| i2(Eo) , 



In the E-M case, the following Sobolev norm for fo = (Eq, Hq): 



So - \\Eo\\ LH v ) + ll V£; o|| L 2 (Eo) + ||V 2 S || L2(So) + ||V 3 ^ || L2(So) 

+ ll#o|| L 4(£ ) + H Vi/ o|| L 2 (So) + ||V 2 i/ || L2(So) + || V3 #o|| L2(So) ■ 

The ideas behind the proof of Theorem [7] are standard. The main points of the 
proof are summarized in further detail in [2lJ Sec. 6.2]. 

Remark. Both Theorem^ and [171 Prop. 6.1] are derived by solving the Einstein 
equations in the CMC gauge along with transported coordinate systems. This, how- 
ever, is by no means an optimal result. For a local well-posedness result requiring 
less differentiability in the CMC gauge in the vacuum setting, see [2] • 

3.4. Outline of the Proof of Theorem [6j Like the other breakdown results dis- 
cussed in the introduction, the proof of Theorem [6] is at its highest level intimately 
tied to the corresponding local well-posedness result of Theorem [7] Our strategy 
for proving Theorem [6] will be analogous to that of the model breakdown problem 
for the nonlinear wave equation |T]). Indeed, we aim to bound the parameters on 
the timeslices in our spacetime which control the time of existence in Theorem [7] 
Assume now the hypotheses of Theorem [6] For each to < t < t\, we define: 

(44) £(t) = ||fc|| L4(ST) + ||Vfc|| £a(Er) + ||V 2 fc|| i2(ST) + ||V 3 fc|| i2(ST) , 

*(t) = \\n L ^ T ) + I|v^|| l2(Et) + ||v 2 ^|| l2(St) , 

as well as a corresponding matter field energy f(r). In the E-S case, we define 



(45) f (r) = HWII^) + ||v 2 ^|| l2(Et) + ||v 3 ^|| l2(Et) + ||vV Hi ^ , 

|V(A<5 



L 4 (S T ) T || v v ^ t ^ll L 2(S T ) 
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L 2 (£ T ) 



+ || V 3 {Ct<t>) 



£ 2 (E T ) ' 



(46) f(r) = ||B|| £ * (Er) + ||Vf;| 




i 2 (S T ) 



+ ||v 3 £ 



iJ IIZ,2(E r ) 

+ llv 3 #||, 



where E,H G X*(M) is the electromagnetic decomposition of F. 
For convenience, we also fix an "initial time" 

max(£o,2fi) < tq < tt, 

and we treat S To as the "initial timeslice" in M. From now on, we will only be 
concerned with the timespan [ro,ti) in this proof. For convenience, we also define 



The mean curvatures of the S r 's, tq < r < t\, are trivially comparable to |ti|. 
Suppose we can also uniformly control the following: 

• The diameters and the injectivty radii of the S T 's, To < t < t\. 

• The quantities &(t), SH(t), and f(r) for all tq <T<t\. 

Then, applying Theorem [7] to each of the timeslices of M, we obtain roughly that 
there exists some sufficiently small < e < |ti| such that for every r < r < t\: 

• A solution to the Cauchy problem exists with initial data given by E T . 

• The solution exists on a time interval including [r, r + e). 

• The solution is unique up to isometric imbedding. 

By combining the above solutions, we obtain the desired continuation of [M, 
to the time interval (to, t\ + e), which completes the proof of Theorem [S] 

As a result, the main objectives will be that of controlling the diameters, the 
injectivity radii, and the energy quantities ^(t), £H(t), f(r). The remainder of this 
paper will be dedicated to establishing these estimates. Of these, the energy bounds 
will be the primary task. The diameter and injectivity radius bounds, on the other 
hand, will be trivial consequences by the end of the proof. 

Rather than directly controlling the quantities associated with &(t), SH(t), and 
f(r), which are horizontal fields on the timeslices, we instead control norms of 
corresponding spacetime quantities on these timeslices. These include the spacetime 
Riemann curvature R and the spacetime matter field, i.e., either (j) or F. The 
necessary estimates for R(t), £H(t), and f(r) can then be derived using standard 
elliptic estimates. For example, R can be related to 1Z and k using the Gauss 
equations (|32[) and the Codazzi equations (|33|) . 

To control R and the matter field we apply variations of standard energy- 
momentum tensor techniques. In particular, we take advantage of the observations 
that both R and $ satisfy covariant tensorial wave and Maxwell-type equations. We 
then construct "generalized energy-momentum tensors" based on these relations, 
and we apply these in the standard fanshion to derive basic energy inequalities. 

The above suffices for "lower-order" a priori energy estimates. In order to derive 
"higher-order" energy estimates, however, we will also need uniform bounds for R 
and for quantities derived from $. Since these quantities satisfy a system of 
covariant tensorial wave equations, we can then apply the representation formula 
presented in Thm. 7] (and also in [HI Thm. 5.1]). 



(47) 



M+ = S [TOitl) = {q e M | r < t (q) < h} . 



; More specifically, D 2 <j) in the E-S case, and DF in the E-M case. 
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The representation formula mentioned above is only valid on the "regular" por- 
tion of past null cones, i.e., prior to the null injectivity radius. Thus, in order to 
obtain satisfactory estimates using this formula, we must also control the geom- 
etry of past null cones. In particular, we must control the null injectivity radius 
and various connection quantities on these cones by other quantities which can be 
controlled a priori. In our case, these include L 2 "flux" quantities for R and $ on 
these cones, and other a priori bounds relating to the time foliation. 

Unfortunately, this task is exceedingly difficult and is responsible for a vast 
portion of the technical work behind Theorem [5] For instance, this involves con- 
structing a geometric tensorial Littlewood-Paley theory, cf. [13] and [2TJ Sec. 2.2], 
and applying it in a massive bootstrap argument. Within this argument are various 
Besov estimates as well as an elaborate sharp trace estimate for regular null cones. 
As a result, we omit a majority of this development from this paper. For details 
regarding this portion of the argument, see [2T]. [3 

Using the above control on the local null geometry along with the representation 
formula of 22! , we can derive all the necessary estimates for R and As a result, 
we can then estimate &(t), 5R(t), and f(r), as mentioned before, and hence Theorem 
|6]is proved. In the remaining sections of the paper, we will provide more detailed 
discussions on the steps described in this outline. 

4. Regular Past Null Cones 

As mentioned before, both the local energy estimates in this paper and the 
representation formula for covariant tensor wave equations in [35] are essential 
components of the proof of Theorem[6l Moreover, both depend heavily on the local 
null geometry of (M, g). More specifically, both are applicable to our setting only on 
"regular" past null cones, where the null exponential map remains a diffeomorphism 
and its image retains a smooth structure. 

We shall provide in this section some preliminaries on such regular past null 
cones. We mainly follow the development given in [3TJ Ch. 3] and [21], but we 
restrict ourselves in this paper to the special case of time foliated null cones. 
This eliminates much of the minor technical irritations present in |21| 122] that arose 
from dealing with more general foliating functions. 

4.1. Regular Past Null Cones. Assume the spacetime (M,g, $), along with the 
time foliation of M given by the time function t, as expressed in the statement of 
Theorem HI Fix p e M, and consider the null exponential map expp about p, i.e., 
the exponential map about p restricted to the past null cone 9t of the tangent space 
T p M. Define the past null cone N~(p) of p to be the image of exp p . 

By definition, N~(p) is ruled by the past inextendible null geodesies beginning at 
p. Recall that exp p is a diffeomorphism between a sufficiently small neighborhood 
of in yt and its image, which is then a smooth null hypersurface of M. We refer 
to such regions of N~(p) as "regular". A loss of such regularity can occur at a 
terminal point z of N~(p), where one of the following scenarios hold: 

• The point z is a cut locus point, that is, distinct past null geodesies from p 
intersect at z. In other words, the map exp p fails to be one-to-one at z. 

1 ^Earlier work in this area for the vacuum case were done in |12l 1131 fLU 1161 [T9l 1241 [25] , 
14 This case was discussed in |2 1 1 Ch. 4] . 

^For convenience, we assume OT does not include the origin of T p M. 
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• The pair p and z are past null conjugate points. In other words, the map 
expp fails to be nonsingular at z. 



Define the function t p on N (p) by t p (q) — t(p) —t(q), i.e., the difference in time 
between the vertex p and the given point q. We also define the following: 

• Let s(p) be the infimum of all values t p (q) for which q G N~(jp) is a past 
null conjugate point. We call s(p) the past null conjugacy radius of p (with 
respect to the t-foliation). In the case that N~(p) has no conjugate points, 
we define s(p) = t(p) — to. Note that s(p) indicates the largest v > such 
that N~(p) encounters no conjugate points before i p -value v. 

• Let \(p) be the infimum of all values t p (q) for which q G N~{p) is a cut 
locus point. On the other hand, if N~(p) has no cut locus points, we define 
\.(p) = t(p) — to- Then, l(j>) indicates the largest v > such that N~(p) 
encounters no cut locus points before i p -value v. 

• Define i(p) = min(s(p), l(p)). We call this the past null injectivity radius of 
p (with respect to the t- foliation). 

We will refer to the region 



morphism between a neighborhood in 91 and J\f~ (p) , and Af~ (p) is indeed a smooth 
null hypersurface of M . From now on, we will only refer to the regular null cone 
J\f~(p), as N~(p) is in general too irregular for our use. 

In addition, for any < v < i(p), we define the null cone segment 



4.2. Normalization and Foliation. Since tangent null vectors in Af~(p) have 
vanishing Lorentzian "length" and are orthogonal to Af~(p), they cannot be nor- 
malized without introducing vectors transversal to Af~(p). Consequently, in our 
treatment, we will require an additional choice of a future timelike unit vector 
t G T P M at p. For this, we use the most natural choice based on the problem at 
hand: the value of the future unit normal T to the E T 's at p. 

We define the null generators of N~(p) (or of N~(p)) to be the inextendible 
past null geodesies 7 on M which satisfy 7(0) = p and g(l'{0), T\ p ) = 1. Wc 
can smoothly parametrize these generators by S 2 using the following process. If 
we choose an orthonormal basis eo,...,e3 of T p M, with eo = T\ p , then we can 
identify each w G S 2 with the null generator 7 W satisfying 7^,(0) = — e +uj k ek- For 
convenience, we assume such a parametrization of the null generators of N~(p), 
and we denote by 7^ the null generator corresponding to u). 

Remark. The objects on M~ (p) that we will discuss are of course defined indepen- 
dently of any parametrization of the null generators. However, for ease of notation, 
we will work explicitly with § 2 . 

In addition, we define L G £(W~(p)) to be the tangent vector fields of the null 
generators of ftf~(p), i.e., we define L\ Ju> ( v ) = l' M { v ) f° r an Y w £ § 2 and < v < x(p). 
We note in particular that L is a geodesic vector field. 

As in previous works, e.g., [121 HH1 EU E21 [24], we would like to express Af~{p) 
as a foliation of spherical cross-sections. This was done abstractly in (STJ Sec. 3.1] 
and (22[ Sec. 2] for the sake of the generalized representation formula for wave 




M (p;v) = {q G Af (p) I t p (q) < v} , 
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equations. Here, we can simplify our presentation by considering only the special 
case of the foliating function t p , the most natural choice for the current setting. 
Define the null lapse function & by the formula 

(48) = (Ltp)- 1 G C°° (AT (p)) . 

Note that $ is everywhere strictly positive, and satisfies the initial limits 

(49) limtfl (v) =n(p). 

The positivity of •& implies dt p is nonvanishing, so the level sets of t p , denoted 

S v = Af~ (p) | t p (q) = v} , v>0, 

form a family of hypersurfaces of Af~ (p) . Since L represents the unique null di- 
rection tangent to AT~(p), and the positivity of i? implies L is transverse to each 
S v , we can conclude that each S v is spacelike, i.e., Riemannian. Furthermore, the 
definition of i(p) implies S v is diffeomorphic to § 2 for every < v < i(p). 

We adopt the conventions of Section |2~21 to discuss this foliation of Af~(p), with 
N = N~(p) and f — t p . In order to distinguish this from the time foliation 
of M, we denote the induced horizontal metrics and connections on the S v 's by 
A and f, respectively. Mixed connections on W~(p) are denoted by y7, and the 
Gauss curvatures of the S v 's are denoted by /C G C°° (J\f~ (p)) . In addition, the 
"normalized normal" vector field in this case is given by Z — i3L, i.e., y tp — ^ l- 

Since J\f~ (p) is null, we have no volume form on J\f~ (p) with respect to which we 
can integrate scalar functions. However, we can still provide a canonical definition 
for integrals of functions over Af~(p). Indeed, we define this integral by 

(50) / (f>= f ( f tf-Adv 

JH-(p) Jo \JSv J 

for any <f> G C co (N~ (p)) for which the right-hand side is well-defined. We can 
similarly define integrals over any open subset of J\f~(p), in particular for J\f~ (p; v). 

Remark. We can show using the change of variables formula that (|50|) is in fact 
independent of the foliating function ofAf~(j>); see \21\ Prop. 3.4]. Later, we will 
further justify (|50p in terms of general local energy estimates. 

4.3. Parametrizations and Null Frames. We can parametrize AT~ (p) using t p 
and a spherical value. For any < v < i(p) and w £ § 2 , we can identify the pair 
(v,uj) with the unique point q on both the null generator 7^ and S v . As a result, 
we can naturally treat any <fr G C°° (J\f~ (p)) as a smooth function on the cylinder 
(0,i(p)) x § 2 . For any such cf>, we denote by <f>\( v . u ) the value of <fi at the point 
q corresponding to the parameters (f,w). Wc will freely use this (v, w)-notation 
throughout future sections without further elaboration. 

In general, null frames are local frames I, rh, ei, ei which satisfy 

9 (j, tj = 9 (m, m) = 0, g (f, rnj = -2, 

g (j, e^j = g (fh, e a ) = 0, g (e a , e b ) = S ab , 

Here, we construct null frames which are adapted to our i p -foliation of J\f~(p). 

Each point of iS„ is normal to exactly two null directions, one of which is rep- 
resented by L. We define L G 3L(Af~(p)), called the conjugate null vector field, to 
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be the vector field in the other normal null direction, subject to the normalization 
g(L,L) = —2. A direct calculation yields the following explicit formula for L: 

(51) L = -n~ 2 -d 2 L — 2n~ 1 i9T. 



Next, we append to L and L_ a local orthonormal frame ei,e2 on the S v 's. Then, 
{L,L, ei,e2} defines a natural null frame for N~{p). We adopt the following in- 
dexing conventions for adapted null frames. 

• Horizontal indices 1,2 correspond to the directions e\ and ei- 

• L corresponds to the index 3, while L corresponds to the index 4. 
Next, if we define the vector field 

N = i (n -1 0L - nd^L) 6 X (Af~ (p)) , 

then (p)Tj) implies g(N,T) = and g(N,N) = 1. In other words, N is the outer 
unit normal to each S v in S t (p)_„. Using (|5ip again, we can derive the following: 

(52) L = -mT 1 (T - N) , L = -n~ 1 d (T + N) . 
With this, we can compute the (/i-) norms of L and L, 

(53) \L\ 2 = 2n a -2 , |L| 2 = 2n" 2 tf 2 . 

4.4. Ricci Coefficients. We will make use of the following connection quantities: 

• Define the null second fundamental forms \, x G TT 2 Af~ (p) by 

X (X,Y)=g(D x L,Y), x(X,Y) = g {D X L,Y) , X, Y e 1 (p)) . 

Both x an d X are symmetric, since both L and L are normal to the Su's. 
We often decompose x into its trace and traceless parts: 

tr X = ^ ab Xab, X = X ~ \ ( tr X) A- 

We also use an analogous decomposition for %. 

• Define C,r? e rrW _ (p) by 

C(X) = i fl (DxL,L), r?(X) = i ff (X,D i L) ! Iel(r(p)). 
From [12, Prop. 2.7], we have the following relation between £ and ij: 

(54) !Z =-C+y(log^). 

The quantities trx, x, and £ are called the expansion, shear, and torsion of 7V~(p). 
We refer to the collection of fields tr%, Xj trx, X, C V as the fficci coefficients of 
A/ -- (p) (with respect to the fp-foliation) . 

From (|51[) and direct calculations, we see that x and are intimately tied to the 
time foliation of M via the following formulas: 



(55) T± a = k la N l +f a {\ogn) , 

(56) Xab = -n- 2 d 2 ■ X ab + 2rT 1 tf ■ k ab . 
Lastly, we define the mass aspect function fj, G C°° (J\f~ (p)) by 

(57) ,i = TCa - \r b X ab + ICI 2 + ~#4343 - \Rm- 



This quantity is present in the representation formula for wave equations, and it 
plays a crucial role in controlling the local past null geometry. 
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4.5. Initial Values. We now briefly examine the initial values of various quantities 
on A/" - (p) , that is, we look at the limits of these quantities at p along null generators. 
This is of importance in the problem of controlling the local null geometry, since 
our goal will be to show that these values differ little from the initial value in 
some specific sense. Moreover, these initial value computations were applied in the 
derivations of the representation formulas of [15j [22] . 

We opt in this paper to skip the technical details on this topic, since most of 
them are rather distant from the heart of Theorem |5J For instance, some of the 
derivations require applications of convex geometry and bitcnsor fields. A complete 
exposition in the case of general foliating functions can be found in |21[ Sec. 3.3]; 
the case of the i p -foliation is covered in [21] Sec. 4.2]. An earlier account for the 
geodesic foliation case, which is needed for the general case, is presented in }24j . 

The first task is to examine the horizontal metrics A at the initial limit. Note 
that a coordinate system (U, ip) in S 2 generates a transported coordinate system on 
each S v by mapping the point with parameters (v, w) to f(oj). 

Proposition 8. Let An denote the Euclidean metric on S 2 , and fix a coordinate 
system (U,(p) on S 2 . Index An using ip- coordinates, and index X on each S v using 
the associated transported coordinate system. Then, for any uj G U , 



lim v Xn\, % = n \ ■ (An) 



l<i,i<2. 



Proof. The proof involves relating the "(i>, w)-parametrization" of J\f (p) with nor- 
mal coordinates; see [21] Prop. 3.15]. □ 

Proposition [8] implies the following integral limit: 

Corollary 9. Let 4> 6 C°°(Af- ip)), let fa € C°°(§ 2 ), and suppose 

<%, w ) = M , we§ 2 . 

Then, the following integral limit holds: 

lim v~ 2 / <ft = n 



,21 



The next proposition deals with the effects of derivatives on initial limits: 
Proposition 10. Let A E FTJ\f~(p), and suppose 

I™ \A\\ {V , U) =0, ^e§ 2 . 

Then, for any integer k > 0, we also have 

\imv k \f k A\\. ,=0, ueS 2 . 

Proof. See j2TJ Prop. 3.17]. □ 

Finally, we consider the Ricci coefficients: 
Proposition 11. The following limits hold for each w £ § 2 . 



• We have the following limits for \ : 
lim Itf(trx) — 2*- 1 1 1 . =( 



(58) lim \ti (tr X ) - 2t^\\ {v uj) = 0, lim |*|| (0|W) = 0. 
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• We have the following limits for £ and r\: 

Proof. Equation ([58]) is an immediate consequence of (2TJ Prop. 3.18]. The limit 
for r\ in (|59j) is an immediate consequence of fj55|) . while the derivation for Q in (|59j) 
can be found in [3TJ Prop. 4.9]. □ 

Remark. We can also obtain analogous initial value properties for both \ and /i. 
For details, see [21] Sec. 3.3, Sec. 4.2]. 

5. A Priori Estimates 

In this section, we derive a priori estimates which are consequences of the break- 
down criterion (|40|) . These are fundamental to the proof of Theorem HI since they 
determine the nature of the null geometry estimates and the higher-order estimates 
which we must prove. In other words, we must control the null geometry and the 
higher-order energy norms by quantities which can be controlled a priori. 

We now introduce the notion of "fundamental constants", i.e., values on which 
all of our "universal" constants will depend. More explicitly, whenever we write 
A < B, we mean A < CB for some constant C depending only on these fundamental 
constants. Similarly, if we write A ~ B, then we mean C~ X B < A < CB for some 
constant C depending only on these fundamental constants. 

The complete list of such fundamental constants is given below: 

• The "breakdown time" t%. 

• The "breakdown criterion" constant Co in (|40|) . 

• The intrinsic and extrinsic geometries of the "initial timeslice" E ro , and the 
values of the fields defined on S To . These include derivatives of R, $, k, n, 
etc., restricted to S To , as well as the volume V(£ To ) of r . 

Remark. We can be more explicit on the exact properties of E To and values on 
S To for this dependence; see |21[ Sec. 6.3]. However, the main idea is that we can 
uniformly control a sufficient amount of quantities ( as dictated in Section \3.J$ on 
all the S T 's, tq < r < t\, by properties of only E To . 

5.1. Regularity of the Time Foliation. Our first task will be to establish some 
basic control for our time foliation. From the breakdown criterion (|4"0)) along with 
([T9]). (|2"T]l. (|23|). (j38)) . we can immediately obtain the following: 



(60) h\\ L oo {M) < 1, HQIU(m) + l|Bic|| £O0(A0 < 1. 

Here, it is the deformation tensor of T, while Q and Ric denote the energy- 
momentum tensor of the matter field and the Ricci curvature of (M,g). 

Contracting (|35[) and recalling the CMC gauge condition yields the lapse equation 

(61) An = n \\k\ 2 + Ric (T, T)l - 1. 

Moreover, the strong energy condition implies that Ric(T, T) > in both the E-S 
and E-M settings. As a result, we can derive the following bounds: 

Proposition 12. On each E r , tq < r < t\, we have the comparisons 

(62) H~l, V(E T )~1. 
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Proof. At a minimum point p of n on £ T , we have from (|61|) that 

i2~ 



Ric (T,T) + |fc| 



> 1. 



By (|60|) . we obtain a lower bound for n depending only on the fundamental con- 
stants. Next, at a maximum point pofnon E T , we have from (|61j) that 



n(p) < 



Ric (T, T) 



y 



where we also applied the strong energy condition on Ric(T, T). 
For a volume form V T on S T , we have 

CtV T = -n (tr k) V T = -nt ■ V T > 0. 

By the above explicit bound n < 3i~ 2 , then 

< -rVCS T ) = -t / re < -3T-V (r) . 
From Gronwall's inequality, we obtain V(r) ~ r~ 3 ~ 1, as desired. 



□ 



Furthermore, from (|40l) and Proposition [T2l we have proved the uniform bounds 



(63) 



li°o(Af + ) + || n |li°°(Af + ) + ll^ 7rl lli°°(A/ + ) + II^IIl°°(M+) ~ 1) 

where A/ + is as defined in (|47l) . 

Next, we examine some coordinate regularity properties satisfied by the times- 
lices. To begin with, since E To is compact, we can find a constant C > 1 and a 
finite covering of S ro by local coordinate systems {U\, ipi), ■ ■ ■ , (U m , (p m ) such that 
for each 1 < i < m, we have the uniform ellipticity bound 

(64) cr 1 ^ 2 < <C|C| 2 , P^.^el 3 , 

where 7 is indexed with respect to the (^-coordinates. We next show that this 
uniform ellipticity condition holds uniformly for all £ T 's, tq < r < t\. 

Proposition 13. Let (Ui,ipi), 1 < i < m, be as given above, and, for each t, we 
let (U[,ipJ) denote the coordinate system on S r obtained by transporting {Ui,<Pi). 
Then, there exists a constant C\ > 1, depending on the fundamental constants, such 
that for any To < r < t\ and 1 < i < m, we have the estimates 

(65) C b - 1 |C| a < 7j -il p ^<Ci,ie| 2 , P eU[,teR 3 , 

where 7 is indexed with respect to the tpj -coordinates. In other words, the uniform 
ellipticity condition (|64p holds uniformly on all the £ r 's. 

We sketch the proof of Proposition [TBI below. For further details, we refer the 
reader to [2TJ Prop. 4.6], as well as [16l Prop. 4.1] and [T7] Prop. 2.4]. 

Proof. Fix a coordinate system (U, ip) — (Ui, tfi) of £ To , and consider its transported 
systems (U T ,(p T ) = (U[,ifl); we will index with respect to these coordinates. Fix 
£ G M 3 , and let X = £%. By (34}, we have C t \X\ 2 = -In ■ k(X,X), so 



\X\ 2 <C t \X\ 2 <2\\n\ 



\X\ 



(66) - 2 Hnll^oo^) II*IL«>(e t ) KM -^i y M ^ ii"-iii~(so ll n 'lli»(2 T , 

Integrating ((66)) along the integral curves of Z — nT (or equivalently, of T) and 
applying Gronwall's inequality, we obtain for each p £U the comparisons 

t el, 



C'- 1 \x\ 2 


< \x\ 2 


< c \x\ 2 




p 


P T 
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where p T is the normal transport of p along Z to E r , and where 



C = exp 



2(«i-t ) \\n\\ Loo(M+) ||vr|| ioo(M+) 



< 1. 



Defining C\, = CC, where C is as in (l64l) . varying over all £ e R 3 , 1 < i < m, and 
p G S To , and recalling the bound (|64|) . then we obtain ([65]), as desired. □ 

By considering the family of transported coordinate systems from Proposition 
fT3l we can deduce some rudimentary Riemannian geometric estimates on the S r 's. 
For a point set 3 and r > 0, we let B(x, r) denote the Euclidean ball in R 3 about 
x of radius r. The property given in the subsequent proposition is an immediate 
consequence of Proposition IT31 and the compactness of the E r 's. 

Proposition 14. Let (U[,ipJ) be defined as in Proposition \13\ for all tq < t < t\ 

and 1 < i < m. There exists ro > 0, depending on the fundamental constants, such 
that for any tq < t < t\ and pe E T , there exists 1 < i < m such that p S U[ , and 

Proof. See [H Lemma 2.2] and [21 , Prop. 4.6]. □ 

Next, for any tq < r < t\, p € S T , and p > 0, we let B T (p, p) denote the geodesic 
ball in S r of radius p about p. In addition, for r and p as above, we define 

. V(B T (p,a)) 

i.e., the volume radius at p (in S T ) with scale 1. By considering the transported 
coordinate systems in the statement of Proposition [131 we can derive the following 
uniform lower bound for the above volume radii. 

Proposition 15. There exists p > 0, depending on the fundamental constants, 
such that r T (p) > p for every tq < r < t\ and peE T . 

Proof. See the proof of [2T] Cor. 2.1]. The main idea is that the uniform ellipticity 
condition implies that geodesic balls are comparable to Euclidean balls. □ 

5.2. Sobolev Inequalities. Proposition [T3l also plays an instrumental role in de- 
riving uniform Sobolev estimates on the E T 's. Indeed, we can derive such first-order 
Sobolev inequalities on each E T by using a partition of unity argument and applying 
the corresponding Euclidean Sobolev inequality to each coordinate system (U[ , tpj). 
As a result of this, we obtain the following scalar Sobolev inequalities: 

Lemma 16. The following hold for any f £ C°°(M), tq <t < t\, and q > 3: 

(67) ll/y (Sr) <HV/|| il(Er) + \\f\\ Ll ^)> 

(68) II/L» (Et) < llV/Hi,^) + ||/|| i4(Er) . 
In (168[) . the constant of the inequality depends also on q. 

Proof. See [12 Lemma 2.1]. □ 



Remark. The significance of Lemma \16\ is that the constants of (|6T|) and (I68|) . 
i.e., the "Sobolev constants", are uniformly bounded over all the S r 's, to < r < t\. 
In particular, they are controlled by the fundamental constants. 



^Also essential to the proof of 1681 1 is the "minimal radius" property of Proposition [14] 
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In addition, we can apply Lemma [16] along with the preceding remark in order 
to derive tensorial versions of Sobolev inequalities. 

Proposition 17. Let tq < r < t\. 

• For any2<p<6, q> 3, and * G TTTM, 

(69) \\n L ^) £ II v*il 2(Et) + ||*|| £3(Et) , 

( 7 °) ii*iIl~( St) <II v *II^( St ) + ||v|/|| ^(^)- 

TTie constants of the inequalities depend also on p and q, respectively. 

• In addition, for any \& G VTJ~M and $ E TTM , 

(71) H*llL- (Er ) ^ llv 2 *!!^) + II*II^( Et) . 

(72) ll^ll^(^) < ||^ 2 *|| l2(St) + P*IIl 2(St) + PII^(e t) • 

Proof. The inequalities (|69|) and f|T0[) are direct adaptations of the proof of the first 
part of [211 Prop. 4.7]; see also (2TJ Prop. 2.3], [2U Prop. 2.4], and [17, Cor. 2.7]. 
The main idea is to apply Lemma [TBI to the scalar quantities \^>\ q . q > 1. lj 

Next, ([?Tj) is proved by applying ([70)1 and in succession. We can also obtain 
([72]) by applying ([70)) and ([69)) in the following manner: 



ll*IL~ (Er) ;$ P*IIl4( Et ) + ll*lli*( ET ) 

^II^ 2 *IL 2 (e t ) + II^IIl 2 ( St ) + II*IIl 2(St )- 

Note that for $ G TTM = TTM, we clearly have |V$| < |D$|. □ 

5.3. Generalized Energy-Momentum Tensors. We now describe some general 
methods for obtaining global and local energy bounds for various spacetime quanti- 
ties. Our focus will be on two general classes of "energy-momentum" tensor fields, 
abbreviated as EMTs, based on covariant wave-type and Maxwell-type equations. 

The "wave equation" EMT was previously applied in [8j [T_7] in order to obtain 
higher-order estimates. In particular, in [T7], this was used to bound the L 2 -norms 
of covariant derivatives of the curvature in a vacuum spacetime. We will apply these 
tensor fields with analogous intentions. The "Maxwell" EMT is an adaptation of 
this process to solutions of "Maxwell-type" equations. We will use this tensor field 
to derive a priori energy estimates for the curvature R, thereby avoiding the more 
precise but more computationally complex Bel-Robinson tensor field. 

Throughout this section, we will assume the following fields: 

• Let U,V € TT r M satisfy the covariant tensor wave equation 

(73) D g U = V, 

where D g is the tensorial wave operator D g — g a ^D a p. 

• Let A G TT r+2 M, B G TT r+1 M, and C G TT r+3 M, with A antisymmetric 
in its first two components and satisfing the Maxwell-type equations 

(74) D a A al jj = Bpi, DyA a pi + D a Apyi + DpA ia i — C^ a pi- 

We begin by describing the "wave equation" case as in [17 . Define the EMT 
Q W {U] G TS 2 M and its associated current P W [U] G X*(M) by 

Qw [U] a p = h {D a U, D P U) - \g a ^ v ■ h (D^U, D V U) , 



We norm mixed tensor fields with respect to h and 7. 
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P W [U] {X) = Q W [U] (X,T), 

where h(-, •) denotes a full /i-metric contraction of the unindexed components. Note 
the resemblance to the energy-momentum tensor (|2"TT) for the scalar field. In ad- 
dition, we define the corresponding divergence quantities D ■ Q W [U] S 3£*(M) and 
D ■ P W [U] e C°°(M) in the expected manner: 

D ■ Q w [U] a = D Q W [U] afj , D ■ P w [U] = D a P w [U] a . 

The main properties of this wave EMT are listed below. 

Proposition 18. The EMT Q W [U] satisfies the following properties: 

• The following estimates hold: 

(75) \Q W [U]\<\DU\ 2 , 

(76) \D ■ Q w [U]\ < \DU\ {\V\ + \R [U}\) + \n\ \DU\ 2 , 
\D ■ P w [U}\ < \DU\ {\V\ + \R [17] |) + \n\ \DU\ 2 . 

• For future causal X, Y £ X(M), 

(77) Q w [U]{X,Y)>0, 
In particular, 

(78) Q W {U](T,T) = ^\DU\ 2 . 

Proof. These follow immediately from a series of direct computations. For further 
details, consult [HI Prop. 4.13] and [17]. □ 

The "Maxwell-type" EMT is constructed analogously. Indeed, we define the 
EMT Q m [A] £ TS 2 M and its associated current P m [A] € X*(M) by 

Qm [A] ap = h (Aa^Ap") - \g af >h (A^,A^) , 
P m [A] (X) = Q m [A] (X,T), 
and define D ■ Q m [A] £ X*{M), D ■ P m [A] e C°°{M) by 

D ■ Q m [A] a = D?Q rn [A] afj , D ■ P m [A] = D a P m [A] a . 

Proposition 19. The EMT Q m [A] satisfies the following properties: 

• The following estimates hold: 

(79) \Qm[A]\<\A\ 2 , 

(80) ^•Q„ l ^]|<| J 4|(| J B| + |C|) + |^||A| 2 , 
\D-P m [A]\<\A\(\B\ + \C\) + \n\\A\ 2 . 

• For future causal X, Y £ X(M), 

(81) Q m [A](X,Y)>0. 
In particular, 

(82) Q m {A](T,T) = ±\A\ 2 . 

Proof. The proofs are analogous to those of Proposition[T8l see |21i Prop. 4.14]. □ 
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5.4. Global Energy Estimates. Next, we apply these generalized EMTs in order 
to establish both global and local energy estimates. By "global energy estimates" , 
we mean L 2 -estimates on entire timeslices, while by "local energy estimates" , we 
refer to L 2 -estimates on past null cones, i.e., "flux estimates". 

We begin with general global estimates. Recall that Killing vector fields can 
be associated via Noether's theorem with conservation laws for quantities derived 
from (standard) energy-momentum tensors. Since we have no Killing fields in our 
setting, we have no hope for achieving such energy conservation. We can, however, 
think of T as an "almost Killing" vector field, since its deformation tensor 7r is 
uniformly bounded. In fact, by slightly modifying the process behind Noether's 
theorem, we can derive energy inequalities by taking advantage of this "almost 
Killing" condition for T along with the generalized EMTs Q w [U] and Q m [A] . 

Lemma 20. Assume U,V, A, B,C as before. If r < T\ < r<z < t\, then 

(83) \\DU\\ 2 l2{ ^ 2) <\\DU\\ 2 l2 ^ ti) + jH jf (\DU\ 2 + \V\ 2 + \R[U]\ 2 )dr, 

(84) UWl^ < \\A\\ 2 l2(Sti) + jH jf (L4| 2 + \B\ 2 + \C\ 2 ) dr. 

Proof. First, for (|83j) . we integrate D ■ P W [U] over the slab Sr TliT2 i. E Applying the 
divergence theorem along with Proposition [5l we obtain the identity 

f Q lv [U] (T, T) — [ Q w [U] [T,T)= [* f n ■ (D ■ P w [U])\ dr. 
By (J63J), G9, dJSJ, then 



\DUWUfr s < \\DU 



2 



L^(S T2 ) ^ H^llL=(s n ) 



1"2 



\DU\ {\V\ + \R[U}\) + \DU 



dr 



Z \\ DU \\l^ rx ) + fj / s (l^l 2 + l^l 2 + \R WW 2 ) dr. 
For ([54^1 . we repeat the above steps with D ■ P m [A\ in the place of D ■ P W [U]. □ 

To obtain the specific a priori estimates we desire, we will apply Lemma [20] to 
the spacetime curvature R and matter field: D 2 (j) in the E-S case, or DF in the 
E-M case. For this, we will need to take advantage of the following facts: 

• The curvature R satisfies Maxwell-type equations. This includes 

(85) D^Rapfos = 0, 

i.e., the Bianchi idenitities, along with the divergence relations (|25|) and 
(|26|) . which hold in the E-S and E-M settings, respectively. 

• In the E-S case, D<j) satisfies the following wave equation: 

(86) a g D a <j) = R ■ D a 4>. 

This can be derived by differentiating (|20[) and commuting derivatives. 

• In the E-M case, F satisfies the following wave equation: 

(87) n g F a p — 2R a Xfl pF\n + R x a F/3x — R pF a \. 

This is derived by differentiating the Bianchi identity for F in (f22l) . 



18 Recall that Sr , = {p £ M \ n < t(p) < r 2 }. 
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For every tq < t < t\, we define the global spacetime energy quantities 

£ l (r) = l|i?ll L2(ST) + ||^|| L2(St) , ^ (r) = PII l2(St) + II^II l2(St) 

in the E-S and E-M settings, respectively. 

Proposition 21. In both the E-S and E-M cases, we have the energy inequality 

& (?)<!, T <T<h. 

Proof. We begin with the E-S case. To handle the curvature, we apply with 
A = R, and with B and C given by ((25)) and (|85|) . in order to obtain 

puw.) < pn^( St0 ) + f / (i^i 2 + i^i 2 i^>i 2 ) 

«/ to ^ 

By (|83|) . with U = Z^, along with (f86j) . we obtain the bound 

\\D 2 4l^ T) < \\D 2 4l {Sto) + jT jf ^ (|^ 2 0| 2 + l^l 2 \Rf) dr'. 

Summing the above equations, recalling the uniform bound (1401) for D(f>, and ap- 
plying Gronwall's inequality to the result yields the desired bound. 

The E-M case is derived in a completely analogous manner. We apply (|83|) and 
(|84|) as before, along with J26|), (|85|) . (|87|) . and we obtain the bounds 




Summing and applying Gronwall's inequality yields the desired inequality. □ 

5.5. Local Energy Estimates. Local energy estimates are derived in a manner 
similar to the global estimates of Lemma 1201 The main difference is that we inte- 
grate over the interior of a past null cone rather than over a time slab. 

Fix p 6 M_|_; we normalize and foliate J\f~(p) as before usin g T \„ and t p . Fix 
also a time value tq < t < t(p) such that t(p) — t < min(i(p), 1). o We define 

N T =N-(p;t(p)-T), l T = {qel-( P )\t(q)>T}, S? = J-(p)nS T) 

where I~(p) is the chronological past of p, i.e., the set of points q G M reachable 
from p via past timelike paths. Note that X T is the region bounded by M T and S r , 
while T,p is the part of E r in the interior of N~{p). 

The following lemma provides the main technical calculation for deriving local 
energy inequalities and justifies our previous definition of integration on M~{p). 

Lemma 22. If lo e X*(M), then 

[ D a u a = I u (T) + / w (L) . 

JX T JY? T jm t 

Proof. See (2TJ Lemma 4.1] or [16]. □ 

"In later sections, we will also define and estimate the higher-order analogues £ 2 (t) and £ 3 (r). 
We adopt the abbreviation £ (t) here in order to maintain consistency of notation. 

2 ^The stipulation t(p) — r < 1 is purely a matter of convenience, in order to eliminate the 
dependence of various inequalities on the timespan t(p) — r. 
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We are now prepared to derive the general local estimate. 
Lemma 23. Assume U, V, A, B, C as before. Then, 

-f q w [u](t,l)<\\du\\I HSPt) 

J N't 

+ C P] L (i^i 2 +i y i 2 +™i 2 ) d<T ' 



(89) -/ Q m [A] (T,L) < ||A|| 2 2(S?) 

(|^| 2 + |B| 2 + |q 2 ) da. 

Proof. Integrate D-P W [U] over X T and apply Lemma W2\ and Proposition[5]to obtain 

/ Q w [U] (T, T)+ f Q w [U] (T,L) = [ ( [ n-D- P w [U]\ da. 

By ([75]) . and ([75]) , we obtain ([55)1 . Repeating the above steps, but with P m L4] 
in place of P W [U] and dHDJ>, ((SJ) in place of (JTHJ), (75]), we obtain flSH|). □ 

Remark. #o<e that both —Q W [U](T, L) and — Q m [A](T, L) are nonnegative. 

Next, we define the lower-order flux densities 

P 1 (p, r) = -Q m [R] (T, L) - Q w [Dfl (T, L) , 

p 1 (p, r) = -Q m [R] (T, L) - Q tt [F] (T, L) 
on A/" T in the E-S and E-M cases, respectively, and we define the corresponding flux 

We can now apply Lemma [23l directly to obtain the desired local energy estimates. 
Proposition 24. In both the E-S and E-M cases, we have the local energy inequality 

•7 7l (p,r)<L 

Proof. For the E-S case, we first apply ([59")) with A = i? to obtain 

- Q m [R] (T, L) < + jy (\\R\\ L ^ a) + \\D 2 4 L2{ ^ } ) da < 1, 

where we also used Proposition [^U Similarly, applying (|88p with U = D<fi yields 



/ Q m [D<t>](T,L)<l. 



The E-M case is handled analogously. □ 

Remark. Estimates analogous to Proposition \ 24\ hold for similar flux quantities 
associated with future regular null cones on the time interval [ro,ti). 

As of now, we have uniform control on the integrals over Af T of some qua- 
dratic quantities in R and either D 2 cf> or DF, given explicitly by —Q m [R](T, L), 
—Qw[D(f>](T, L), and — Q W [F](T, L). The final step here is to determine exactly 
which components of R and the matter field are controlled by Proposition [Ml 
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Proposition 25. OnM T , we have the comparisons 
-Q W [U] (T,L)~ |F 4 C/| 2 + |Ff/| 2 ; 

2 2 

-Q m [A](T,L) ~\A 43 \ 2 + ^|^ 4a | 2 + ^ \A ab \ 2 , 

a—l a,b—l 

where we have indexed with respect to adapted null frames on Af~ (p) , and where 
A a p G YT r M is defined {A a g)i = A a pj. 

In other words, —Q W [U](T,L) controls all the components of DU except for 
those of the form D3U1, i.e., derivatives transverse to Af T . Similarly, — Q m [A](T, L) 
controls all the components of A except those of the form A^a] . This coincides with 
standard results for scalar and Maxwell fields in Minkowski space. 

Proof. We begin by defining the quantity tp = g(T, L), which satisfies 

H™, f\(v u0 = 1 > \ C t p f \ = ntp- 1 ■ \Lip\ < nip- 1 \tt\ \L\ 2 < <p, 

where we have applied ([5^|. Integrating the above along the null generators of 
N~{p) and applying Gronwall's inequality yields the comparison tp ~ 1. E3 

Now, the proof is simply a matter of expanding — Q W [U](T, L) and —Q m [A](T, L) 
in terms of adapted null frames. From these computations, we obtain 

-Q w [U] (T, L) = \p- x |f 4 U\ 2 + \p |f f/| 2 , 

-Q m [A}(T,L) = l<p\A 43 \ 2 + ~p- l '£\A 4a \ 2 + ~<p \Aab\ 2 - 

a—l a, 6—1 

For some more details on these computations, see [21] Prop. 4.17]. The desired 
results now follow from the comparison tp ~ 1. □ 

We also note the following simple consequence of the proof of Proposition |25| 
which will be essential in later estimates involving past null cones. 

Corollary 26. On M T , the null lapse -d satisfies 1? ~ 1. 

Proof. Since 1? = ntp^ 1 , where p is defined as in the proof of Proposition 1251 then 
the desired estimate follows from (1551) and the estimate p ~ 1. □ 

Remark. As a result, the quantity $ will be negligible in estimates. In particular, 
the quantity ±U in Provosition \25\ can be replaced by)P tp U. 

Lastly, we use Proposition [25] to see which null frame components of R and 
the matter field are controlled using Proposition |2~41 First, for R, Proposition [231 
implies that the components which are not controlled are those of the form i?3 Qa /3 . 
However, using the symmetry R^ a af3 = Ra/33a, we see in fact that the uncontrolled 
components are those of the form i?3 a 3t>. More explicitly, if tp denotes a null frame 
component of R not of the form i?3 a 3b, then 

IMI^(Af T ) S 1. 

Remark. The same result could also be obtained for R by using instead the Bel- 
Robinson tensor field. This was done in [161 117] in the E-V setting. 



Here, we also take advantage of the assumption t(p) — r < 1. 
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Next, in the E-S case, Proposition 1251 implies that only the components D^ a (j) of 
D 2 4> are not controlled. Since D 2 tfi is symmetric, then in fact, only D^(f) cannot be 
controlled; for any other null frame component ip of D 2 <fi, then 

Finally, in the E-M case, Proposition l25l implies that only the components D^Fa^ 
are not controlled. By rearranging components via the Maxwell equations, we can 
easily establish that the only truly uncontrollable components of DF are those of 
the form D^F^ a ] for any other null frame component tp, then 

Remark. The second fundamental form k also satisfies a tensorial covariant wave 
equation. One can also use this equation to derive analogous global and local energy 
estimates for k. Such global energy estimates are equivalent to those obtained from 
the elliptic estimates of (196|) . while the local estimates are a new result. 

5.6. Preliminaries for Elliptic Estimates. The last class of a priori estimates 
we will establish are elliptic estimates on the E T 's. Before beginning this task, we 
first focus on some general concepts which will be used in these derivations. 

A primary idea behind these elliptic estimates is the relations that exist between 
corresponding horizontal and spacetime objects, e.g., the curvatures 1Z of the E T 's 
and R of M. Our next step is to briefly describe such relations in an informal 
schematic manner. A more detailed account can be found throughout |21j . 

For any A £ TTM, we let T1A £ TTM denote its projection to the £ T 's. In 
particular, if A is a vector field, then UA is its orthogonal projection onto the 
tangent bundles of the S r 's. If A is fully covariant, then 11^4 is given trivially by 

UA(X 1 ,...,X n )=A(X 1 ,...,X n ), X ll ...,X n eX{M). 

Furthermore, if A is scalar, then HA = A. 

For simplicity, we will adopt an informal schematic language to describe tenso- 
rial expressions in an inexact form. We describe the conventions here merely by 
example; for a more exact description of the schematic notation, see |21j . Consider 
the Gauss equations (|32|) . which yields the informal schematic relation 

(90) K « ILR + k ■ k. 

Here, the symbol • represents a tensor product followed by a finite number of 
contractions and metric contractions. Moreover, the multiple terms of the form 
k ■ k in (|32|) are compacted into a single term in the informal equation (|90"]) . 
We take advantage of such schematic relations by using the following lemma: 

Lemma 27. //$ £ TTM is fully covariant, then we have the schematic relation 

(91) VIFF « ILD$ + k ■ n ($ • T) . 

Proof. Suppose $ £ TT°M, X £ X(M), and Z £ YT r M . Then, by definition, 

V x (n$) (Z) = X ($ (Z)) - $ (V X Z) = D X <I> (Z) + $ (D X Z - V X Z) . 

The first term on the right is simply n£>$ applied to X®Z. Moreover, the quantity 
Dx Z — VxZ can be expressed as a sum of terms of the form k ■ X ■ Z ■ T. ri The 
desired equation ((9~ij) follows as a result of the above observations. □ 



; For example, if Z e £(M), then D X Z - V X Z = -g(D x Z,T)T = -k{X,Z)T. 
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For instance, considering (|90|). if we wish to control VIZ, then applying Lemma 
[27] with $ ps R, we need only control terms of types DR, k ■ R • T, and V/c • k. El 
As a result of the above, we obtain the derivative bound 

(92) \Vn\<\DR\ + \R\\k\ + \Vk\\k\. 

With a little more care, we can iterate the above argument using ([90]) in order to 
estimate higher derivatives of 7Z. In fact, such estimates for V 2 1Z will be necessary 
later when we establish higher-order elliptic estimates. 

In addition to Lemma (37] we will a l so make use of the following general elliptic 
and Bochner estimates. For any fully symmetric £ G TTfM, s > 0, we define 

v^err^M, (v-0j = v°&i, 
v x £ e rr +1 M, (v x = Va6i - v 6 £ a /, 

'(tr0 7 = 7^«M *>1, 
v s = l, 

i.e., the divergence, curl, and trace of £, respectively. 

Lemma 28. If £ £ TTM is fully symmetric and nonscalar, then 

(93) ' - ' "~ ' * ^ ~ 



tr£ e rr s ~ 2 M, 



£ ivei 2 < jf (iv.£i 2 + i|vx£i 2 + i^i£r 

Moreover, for any * e TTA/ and £ C°°(M), 

(94) / |V 2 *| 2 </ ("|A*| 2 + |V*| 2 + i |^| 2 1*| 2 



|V 2 </>| 2 < / (|A0| 2 + |ft||V0| 2 

Proof. See [2TJ Lemma 2.3, Cor. 2.5]; see also Ch. 2]. □ 

In particular, Lemma[28]will be useful for elliptic estimates on k and n. Variants 
of Lemma [28] will also be essential for controlling the local null geometry. 

5.7. Elliptic Estimates. Recall that our primary task in proving Theorem[B]is to 
control the horizontal energy quantities R(t), 9t(r), and f(r) defined in (|44]) - (|461) . 
To accomplish this, we must apply various elliptic estimates in conjunction with 
energy estimates derived using EMTs, for example, Propositions [21] and [24] In 
addition, we must control L 2 -norms of derivatives of the lapse n, as well as obtain 
an L°°-bound for Ctn. These will be necessary in order to control the geometry of 
null cones and derive higher-order elliptic estimates. 

We begin here by deriving lower-order a priori elliptic estimates for k and 1Z, 
as well as the horizontal formulation f of the matter field. These estimates are the 
E-S and E-M analogues for the elliptic estimates of |17[ Sec. 8] in the vacuum case. 
The main difference here is that we must also bound the matter field. 

Remark. The proofs of \17\ Sec. 8] were adapted for maximal foliations. Here, we 
simplify the process by taking advantage of the volume bounds of Proposition \1S\ 

Recall that the "horizontal matter field quantities" f are given as follows: 
• In the E-S setting, this consists of the scalar field (j> itself and L t <f>. 



23 Note that |IK>| < |*| for any $ G TTM. 
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• Iii the E-M setting, this consists of the pair E,H € X(M) representing the 
electromagnetic decomposition of F. 

Remark. Recall that given a set orientation of M (or equivalently, ofT, T[) ), the 
electromagnetic decomposition of F is given by the formulas 

(95) E t = F m T a , H t = *F m T a = -e^ F^T<* . 

First, by Proposition [T2l and (j40]). we have the following trivial bounds: 
Lemma 29. For any 1 < p < oo and tq < r < t\, we have the estimates 

\\ k \\ L p(s T ) + IIv</>|| lp(St) + \\^t<P\\ LH s T) < i, 
II*IIl*( Et ) + II^IIlp(s t ) + \\ h Lp(x t ) Z i 

in the E-S and E-M cases, respectively. 

The next batch of elliptic estimates are consequences of various schematic rela- 
tions between horizontal and spacetime objects. 

Proposition 30. For any r < r < t\, we have the elliptic estimates 

(96) \\n L ^ T) + HVfc|| ia(Er) < 1. 
Moreover, for any t < r < t\, the estimates 

(97) ||V 2 0|| l2(St) + ||V(£^)|| l2(St) <1, 

l|VB|| La(Er) + ||VJf || La(Er) < 1 

hold in the E-S and E-M cases, respectively. 
Proof. From (|32[) . we immediately obtain 

II^IIl=(e t) ;$ II^IIl= ( e t ) + l|fe|li»c=r) ~ 1 + H^IIl^o • 

To bound Vfc, we appeal to (j93| . From (|33|) and the CMC gauge condition, we 
obtain |V • k\ + V x k\ < \R\, hence by (|93[) . Lemma [29l and the above, 

l|v*|| £3(Er) < l|i?ll L2(ST) + (/ s l^i |fc| 2 ) 2 < l + \\r\\ LHSt) ■ 

Combining the above and applying Proposition [2T1 yields (|96l) . 

For the matter field estimates of (|97|) . we appeal to the informal schematic 
language. In the E-S setting, we have the schematic relations 

(98) V0«ILD0, n- 1 £ t 0«n(D0-T). 
Similarly, in the E-M setting, 

(99) E^U(F-T), H «n(V -F-T), 

where V is the volume form for M (given an orientation of M). The estimates 
(|i?T)) then follow immediately after applying Lemma |2"T1 and Proposition [5T] to the 
schematic equations (|98|) and (j99|) . For details, see [211 Lemma 6.5]. □ 

We can use the lapse equation (|61j) in order to derive some higher-order a priori 
bounds for n. These are described in the next two propositions: 
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Proposition 31. For any to < r < ti, we have the estimates 
(100) 



C t n 



< 1, 



(101) -AM+ +Ric (T,T) 

(102) l|V(>C t n)|| La(ST) + ||£ t n|| L2(Er) < 1. 

Proof. Equations PD|) . (|61l) . and Proposition IT21 imply ||An||^oo(2 x ) ^ 1- Then, 

/ |V 2 n| 2 </" (|Ati| 2 + |V ? i| 2 ) < 1 + / \K\ 2 , 

by ((911) and Lemma [29l applying ([96]) results in ([TOO]) . 

The proofs of (jlOip and (|102[) are a bit more involved. Recall that £t applied to 
a scalar tfi is the same as nTcf>. We begin by commuting Ct with A to obtain 

\AC t n - C t An\ < \n\ \k\ \V 2 n\ + (|Vn ■ k\ + \n ■ Vfe|) |Vn| , 

where we also used (|34|) . Applying (j40]l . Proposition [30l and (|100j) yields 

(103) ||AAn-AAn|| i3(Sr) <l. 
Next, differentiating (|5T|). we have 



£ t An = 



|fc| +Ric(T,r) £ t n + 7i-£ t |fc|' 



£ t [Ric (T,T)] 



Moreover, applying (|35j) and (|63|) . we obtain 



n • C t \k\ 



< 1 + \R\ 



|V 2 ?i| 



If we let $ denote Dtp in the E-S case and F in the E-M case, then by (|19|) and 
(|40|) . we sec that |DRic | < \D$\. As a result, 

\n-C t [Ric (T, T)]\ < |DRic| + \D T T\ |Ric| < 1 + \D$\ . 

Combining Proposition |2"T1 (|100p . (|103p . and the above, we obtain (jlOip . 
Define now the scalar 



P = -AC t n 



Ric (T, T) 



-tn. 



Multiplying P by Ctn, integrating by parts, and decomposing \k\ 2 , then 



St 



7*1 



< 



iv£ 



Ric (T, T))\C 



-tn\ 



\P\\£tn\, 



where we have also used the strong energy condition satisfied by both the E-S and 
E-M models. Applying a weighted Cauchy inequality to the right-hand side along 
with ([TUT]) yields (fTP2|) . as desired. □ 

Proposition 32. For any tq <T<t\, 

P 2c Hl^ t) + HAn|U (Sr) < 1. 
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Proof. By (j40j). (|94|. and (|96|). 



|V 2 An|| i2(ST) < \\AC t n\\ m 



-AC t n 



li 4 (S T ) 

\k\ 2 + Ric (T, T) 



<I|£H| l2(Et) + 
+ I|V£^|| L 4 (St) . 

Applying Proposition [31] then 

ll V2 ^llL 2 ( ST )- 1 + l|V£tn|l ^(^)- 

Next, by the Sobolev estimate [T71 Cor 2.7], along with (11021) . then 

'I V 2 £ t n|| r2l ^ , < 1 + || V 2 An||? 2 ^ , ||VAn|| J 2(St) + || V£ t n|| L2(s) 



i 2 (S T ) 



li 2 (S T ) 



< 1 



li 2 (E T ) 

■e3 ||V 2 £ t n| 



li 2 (S T ) 

for any e > 0. Taking e sufficiently small, we obtain the desired bound for V 2 £tn. 
Lastly, we can bound Ctn in L°° by the Sobolev bound (fTl]) and the above. □ 

The uniform bound for Ctn in Proposition l32l omitted in the earliest versions of 
the works in the E-V setting, is essential for controlling null injectivity radii. 



6. Local Regularity of Null Cones 

In this section, we will summarize what is the most technically demanding step 
in the proof of Theorem [6) control of the geometry of regular past null cones. This 
step will be crucial for applying the generalized representation formula for tensor 
wave equations, which will be needed in order to obtain the higher-order energy 
bounds required to complete the proof of Theorem [5] 

Throughout this section, we fix an arbitrary point p G M+, and we normalize 
and foliate the regular past null cone Af~(p) by T\ p and t p = t(p) — t, as usual. In 
addition, we fix a constant < <5o < min(i(p), 1), and we define the segment 

N = N- (p; S ) = {z e AT' (p) | t p (z) < 6 } . 
The main result of this section states roughly the following: 

• The past null injectivity radius i(p) olp is bounded below by some constant 
5 > depending only on the fundamental constants. 

• For small enough 80, depending only on the fundamental constants, the 
Ricci coefficients x, X> C> V °f Af~(p) (see Section |4~4|) can be controlled on 
M in various norms by the fundamental constants. 

The explicit result is stated in Theorem [33l and a detailed proof of the theorem is 
available in [2TJ Ch. 7-8] . Due to the prohibitive length and the level of technical 
detail of the argument, we omit the proof of this theorem in this paper and refer 
the reader to [5T] for details. El On the other hand, we do provide a brief outline of 
the main components of the proof at the end of this section. Our focus, however, 
will be on the consequences of Theorem [33] in relation to the breakdown problem. 



See also 1 1 21 1141 1161 1191 1241 125] for earlier results in Einstein- vacuum settings. 
2 ^In fact, a complete accounting of the details of the proof of Theorem 1331 would more than 
double the length of this paper! 
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6.1. Motivation and Past Results. In Riemannian geometry, one generally re- 
quires L°°-bounds on the curvature in order to derive uniform lower bounds for the 
injectivity radius. Similarly, for the Lorentzian case, given an L°°-bound for R, we 
can control the Ricci coefficients x.i Xi Ci an d V without too much effort. Uniform 
bounds for the null injectivity radius would also follow. 

Unfortunately, we will only have the local flux bounds of Proposition [24] here, 
which makes our task tremendously more difficult. Indeed, bounding the Ricci 
coefficients using L 2 rather than L°°-estimates will necessitate the use of sharp 
trace estimates along the generators of null cones, which then requires the geometric 
Littlewood-Paley theory of [13] and the resulting Besov estimates. This process is 
responsible for much of the technical difficulties in the breakdown problem. 

Such i 2 -curvature results were first obtained by S. Klainerman and I. Rodnianski 
in [T^l HH US] . The results of [HJCH], however, applied only to truncated null cones 
and hence were not directly applicable to the breakdown problem. To address this 
issue, Q. Wang, in [24] [25], extended the estimates in [T2] [14] to past null cones 
with vertex initial data by keeping track of scaling factors at every step. Later, 
D. Parlongue, in [19], revisited this argument for time-foliated null cones. 

In all the above, an essential assumption is that the spacetime is vacuum. The 
primary contribution of this section is the extension of this family of results to 
the E-S and E-M settings. We offer now a more detailed comparison between the 
current results and its predecessors: 

• The papers [12l[T4][T9] all considered the case of truncated regular null cones 
with prescribed spherical initial data. In contrast, we adopt the setting of 
[24l [25] and consider the case of regular null cones with initial data given 
by a vertex point. As mentioned before, this is the type of result needed 
by the breakdown problems in this text as well as in [17] [19] . 

• Unlike [T2] [T4] [24] [25] , which bounded the Ricci coefficients with respect to 
the geodesic foliation, we obtain these bounds in terms of the time foliation, 
like in [19]. This has the advantage of being able to interface directly 
with the breakdown problem. In this "gauge" , we will need the regularity 
properties of the time foliation in order to bound the Ricci coefficients. 
However, this also results in improved estimates for both r\ and X- 

• While pH HU QU [HI (24] ES] dealt exclusively with the Einstein- vacuum 
setting, here we extend the results to E-S and E-M spacetimes. In par- 
ticular, the "curvature flux" defined in QU UM EH (Ml 123 is replaced by 
an analogous flux quantity involving both the curvature and the matter 
field; see Proposition [24] Furthermore, several of the structure equations 
governing the Ricci coefficients now contain additional terms reflecting the 
contributions from the nontrivial matter field. 

• The works [HI HU HH [Ml ES] all dealt only with the special case of very 
small curvature flux and initial values, over a unit interval along (possi- 
bly truncated) null cones. This is, however, not directly applicable to the 
breakdown problem, since in this setting we can only stipulate bounded, 
not small, curvature flux. In our case, Propositions |2~T1 and 1241 only estab- 
lished the existence of some possibly large a priori energy and flux bounds. 
On the other hand, we can in the breakdown problem work in only a very 



2 ^The works 1241 1251 also addressed an error in |12| . which assumed a Besov bound that does 
not hold. In fact, a significant amount of effort in 24 25 is dedicated to addressing this issue. 
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small interval along null cones. Heuristically, one can see that this "large 
flux, small interval" case relates to the "small flux, unit interval" case via 
a rescaling argument. However, we handle the "large flux, small interval" 
setting directly in the statement and proof of Theorem 1331 
• The only paper in the existing literature to address the null injectivity 
radius is |16) , which was separated from the remaining components [12[ 114] 
in the overall proof. In the time foliation case, however, the interplay 
between the "Ricci coefficients" and the "null injectivity radius" portions 
of the proof is more subtle and must be addressed in tandem. 



6.2. Integral Norms. In order to state the main result, we must first define the 
relevant integral norms used within. We begin with natural integral norms on M 
(or Af~(p)): for any 1 < q < oo and * 6 rTT7V~(p), we define 



(104) 



1*1 



1*1 



where mixed tensors are normed with respect to h and A, and where the above 
integral over M is defined using the formula (|50p . In addition, wc will make use of 
the following "null Sobolev" norm on J\f: 



(105) 



1*1 



|y*| 



L 2 (J\f) 



\t~ 1 ty\\ 



H 1 (AO - l|y*p*llz,2(AQ 

Given 1 < q < oo and 1 < r < oo, we also define the following iterated norms: 

i 



(106) 



L q t L r (Af) 



\nir {Sv) dv 



1*1 



L?° L r {J\f) — SU P 
0<i><(5 



1*1 



As before, we use S v to denote the level set {z e JV t p (z) — v}, while L r (S v ) 
refers to the tensorial L r -norm on S v with respect to the induced metric A on S v . 
Since ~ 1 on M due to Corollary [Ml then by ([50]). 



1*1 



LlLi(N) 



1*1 



Li(M) 



1 < q < oo, 



i.e., the L\L q and L 9 -norms on M are equivalent. 

Recall that M can be parametrized by the ip-value < v < So and a spherical 
parameter ui £ S 2 . We can then define additional iterated norms with respect to 
this parametrization: for any 1 < q, r < oo, we define 



(107) 



1*1 



So 



dv 



l*lltr££(A0 = sup / l vE 'l r l(-,-) da; 

* " v 0<v<6 \Js, 2 



1*1 



JO \uj£S 2 / 



dv 
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We can reverse the order of integration to obtain additional useful norms: 



(108) 



1*1 




1*1 



I* 



/ ( sup |*|| ( 
JS 2 \0<v<8 



I if (AT) = S U P \J 1*11 dw 



did 




Lastly, the LfL^{N) and L^Lf 3 (AO-norms coincide with the L°° (AO-norm. 

Later, we will see that the parameter L^L^-norms are comparable to rescalings 
of the natural iterated L(L r -norms. Of the reverse parametrized norms, the most 
important will be the "trace" L^L^-norms and L£Lj° -norms. 

6.3. The Main Result. With all of the appropriate norms defined, we can now 
state the main result of this section. 

Theorem 33. For any q S M+, the following hold: 

• There exists a constant So > 0, depending only on the fundamental con- 
stants, such that i(q) > mhi(5o,t(q) — To). 

• Letting J\f q = M~{q; min(5o, t(q) — To)), then the following estimates hold: 

2 



(109) 



0(tr X ) 



t„ 



0(trx) 



2 



+ IIxIIlp 



L 2 t (M q ) 



lien 



< l, 



< i, 



where x and £ refer to the corresponding Ricci coefficients of AT (q). 
In addition, for Af q as above, we have 



(110) 



^ (trx) - - 

tq 



3 

till 



where fi is the mass aspect function on Af q defined in (|57jl . 
Remark. An analogue of Theorem \33\ holds for future null cones N + {p). 

With greater care throughout the proof of Theorem [33l we could in fact be more 
precise about the bounds for the Ricci coefficients in ()109|) and (|110[) . and how they 
are affected by the fundamental constants. Such a task is significantly simpler in 
the geodesically foliated E-V settings examined in [T^l [2H [55] , since in these cases, 
the only external parameter in the problems is the curvature flux. 

Again, the details of the proof of Theorem l33l are left to [2TJ Ch. 7-8] due to its 
length and the amount of technical background involved. A brief summary of the 
proof is given at the end of this section, after first discussing the consequences of 
Theorem l33l pertaining to the current breakdown problem. 
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6.4. Some Basic Consequences. We now list some basic consequences of The- 
orem [33] which will be useful in upcoming analyses involving null cones. Recall 
that the induced metrics on the S v 's are denoted by A. We define the "rescaled" 
metrics on these level sets by A = (t p )~ 2 \. This rescaling is essential, since in the 
case of null cones with vertex initial data, the uniformities occur with respect to A 
rather than A. For example, from Proposition [8j we observe that A tends toward 
the standard Euclidean metric on § 2 (rescaled by a factor of n(p) 2 ) as t p \ 0. 

One consequence of Proposition [8] and Theorem [33] then, is that as long as the 
timespan So of TV is sufficiently small, then A does not differ much from the standard 
Euclidean metric on § 2 . Another related result is that for similarly small So, then 
the volume forms V on the <S„'s with respect to A do not differ much from that of 
the standard Euclidean metric on § 2 . Precise statements, in terms of transported 
coordinate systems, can be found in [211 Lemma 7.1, Lemma 7.3]. 

This latter property immediately implies the following integral comparisons: 

Proposition 34. Suppose So < i(p) is sufficiently small with respect to the funda- 
mental constants. Then, for every < v < Sq and VP £ YT1~N ' , we have 

m~V 2 I \*\\ M db>. 
In particular, for any 1 < q, r < oo, we also have 



tp-9 



LlLUU) > 11*11^^11*11^.(^0 



In other words, the natural iterated L\ ZZ-norms are equivalent to the (v, u)- 
parametrized norms. The latter representation is especially important, since it 
is expressed as an integral over the product (0, Sq) x § 2 , and hence the order of 
integration can be easily reversed. These reversed L^Lj-norms will be essential in 
numerous estimates, especially within the proof of Theorem 1331 

The regularity of the metric A, as described in [21[ Lemma 7.3], can be used to 
derive uniform first and second-order Sobolev estimates in a process analogous to 
that of Proposition [T71 The results are given below: 

Proposition 35. Suppose 5q < \{p) is sufficiently small with respect to the funda- 
mental constants. Then, for every < v < So, 2<r<oo, and ^ £ TTM , 

ii*iiL«^)<« 1 ^iiy*ii^( S „)+«-" 11*11^), 
\\n L r {Sv) < wn^fsv) n*ii!* (5 „) ii*ii 

<^W^\\l 2 (s v )+ v ^ 1 11*1 W„) > 

^ v \\f 2 Hms v) + v ~ 1 11*11^(5.)- 

In the first two estimates, the constants of the inequalities also depend on r. 

Proof. See [HI Sec. 2.1, Sec. 7.2]. □ 

Remark. The powers of v present throughout the estimates of Proposition^^ are 
consequences of the fact that the uniformities of this setting occur with respect to 
the rescaled metrics A rather than the induced metric A. 
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In addition to the Sobolev inequalities of Proposition I35[ which are estimates 
on the individual S v 's, we have the following first-order "null Sobolev" inequality, 
which is in contrast a Sobolev-type estimate on all of J\f. 

Proposition 36. Suppose 5q < is sufficiently small with respect to the funda- 
mental constants. Then, for any * G VT\N ', the following estimate holds: 



1*1 



a* 



< 11*1 



As a consequence, the following estimates hold for any 2 < d < oo: 



Proof. See [3TJ Lemma 7.8]. 



to' * * 



< 11*1 



□ 



In Theorem [33j we have sufficiently controlled the Ricci coefficients \ and C- The 
remaining coefficients x an d V can also be controlled by observing their relationships 
with x, £, and the quantities fc and n derived from the time foliation. 

Proposition 37. Suppose Sq < \{p) is sufficiently small with respect to the funda- 
mental constants. Then, the following bounds hold on At ': 



^(trx) + - 



Proof. The L°°-bound for n follows immediately from (|40| and ([55]). and the trace 
bound for y/-d follows from this due to (|54|) and (|109[) . To handle the estimates for 
X, we resort to the identity (|56[) . from which one can derive the bounds 



■& (tr x) 



< 1 



tf(trx) 



Ixl < 1 + Ixl 



on AT. The desired estimates for x now follow from f!109[) and (|II0l) . 



□ 



Remark. It is also possible to obtain estimates for r\, ^d, and x in the H^-norm 
analogous to those for x ond £. Such estimates are in fact important within the 
proof of Theorem \33\ itself. However, we will not need these results in this paper. 

Remark. If K, is the Gauss curvatures of the S v 's, then the quantity K — ^~ 2 £~ 2 
also satisfies estimates in both the L 2 (Af)-norm and the H~ 1 / 2 {J\f)-norm. The 
latter bound is particularly essential to the proof of Theorem \33\ since it validates 
numerous elliptic estimates on the S v 's that play a fundamental role in the proof; 
see |21[ Sec. 2.3, Sec. 7.2]. Since these bounds play no role outside of the proof of 
Theorem \33l however, we omit serious discussions of this topic within this paper. 

In the remainder of this section, we will sketch the proof of Theorem [33] For the 
sake of clarity, this will be divided into multiple steps. Again, the reader is referred 
to [H] for a detailed account of the proof. 



27 See [21] Cor. 4.3]. 
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6.5. Proof Outline I: Bootstrap Assumptions. The main argument is a "boot- 
strap" , in which we assume much of what we are trying to prove and proceed to 
derive even better estimates. We can then establish via a standard continuity ar- 
gument that these estimates hold even without the assumptions. 

To be more specific, we define the following two conditions, which correspond to 
the conclusions of Theorem l33l Let p e M+, and let < <5o < 1 and Ao > 1 be 
fixed constants, whose values will be determined later. 

• (NO) PySQyAo : If 5 = mm(S ,t(p) -r ,i(p)) and Wo = Af~(p;6), then 



(111) 



An > 



tf(trx)-- 



tf(trx)-- 



+ WMl^L^ATo) + llxll«i(A^ ) + IICII L ~L2 (M)) + ||Cll^i(JV- ) 



mi 



K 



mi 



• (Nl) Pi j 0i A : If 8 = min(<5o,£(p) — To), then \(p) > S, and the inequality 
(jllip holds for the null cone segment Ao = Af~(p; 5). 

Note that the (N1) Pi< 5 0i a condition implies the (NO) Pj( 5 0j Ao condition. Further- 
more, the conclusions of Theorem l33l imply that the (Nl) gj( $ 0: A condition holds 
over all q £ M+ for some < Sq < 1 and Ao > 1 depending only on the fundamen- 
tal constants. The (NO) p ,,5 0i Ao condition, on the other hand, corresponds to the 
conclusions of Theorem 1331 except the null injectivity radius bound. 
We begin the proof of Theorem [33] by assuming the following: 

The conditions 0^O) q j o _ / \ Q holds for every q 6 M + , with the associated 
parameters < <5o < 1 "sufficiently small" and Ao > 1 "sufficiently large" , 
both depending only on the fundamental constants. 

This is the explicit bootstrap assumption we use for our proof. The precise require- 
ments for Sq and A in our bootstrap assumption is determined within the proof. 
Our goal is then to show the following estimates: 

The conditions (N0) q ^ (u A /2 hold for every q G M + . 

Once we establish this, we can conclude that the conditions (N0) g ^ ^ hold even 
without the bootstrap assumptions. El This is essentially the bootstrap argument. 

To carry out this process, we require numerous auxiliary estimates, by which 
we mean estimates that are consequences of the (N0) Pi 5 .a condition. This is 
fundamental to the proof, as many of these auxiliary estimates imply sufficient 
regularity on the S v 's so that the relevant analysis tools can be applied. 

The primary examples of auxiliary estimates are the properties described in 
Section WM Note that these properties, in particular the bounds and comparisons 
of Propositions l3~ill3"Tl depend only on the (N0) p ,5 0i a condition and not on the 
validity of Thcorcm[33litself. For example, the Sobolev estimates of Propositions [35l 
and[36]are vitally important for various basic geometric Littlewood-Paley estimates. 

Another important class of auxiliary estimates are the symmetric Hodge-elliptic 
estimates on the S v J s. These arguments are essentially variations of those originally 
introduced in [71 Ch. 2]. First, we define the following operators: 



28 We of course also need the initial value results of Proposition 1111 
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• Given a horizontal 1-form £ 6 X(Af) , we define 

e C°° (Af) x C°° (AO , T> X Z = (\ ab f a Z b , e ab f a £ b ) , 

where e is the associated volume forms on the S^'s. 

• Given a symmetric traceless £ £ T 2 Af, we define 

• Given G C°°{N) x C°°(Af), we define 

The notation T>\ is justified by the fact that D\ is the L 2 -adjoint of T>\. 

Here is where the L?°H _1 / 2 -estimate for K,—d~ 2 t~ 2 becomes essential. From this 

t 

assumption, we can prove that T>\ and 2?2 are one-to-one, and hence we can define 
via projections their L 2 -bounded "inverses" Vf 1 and T>2 1 . Although T>\ clearly 
fails to be injective, we can still define a viable "inverse" (P^)" 1 by "modding out" 
constant functions and mapping to the canonical element with zero mean. For 
detailed explanations of these procedures, consult [HJ Sec. 2.3]. The upshot of this 
development is the following set of auxiliary Hodge-elliptic estimates: 

Lemma 38. If the (T^O) Pt s .A condition holds for sufficiently small Sq (with respect 
to Ao and the fundamental constants), then the following estimates hold: 

II^^L=(5„) 11^11^(5.) S MWlHS.) > V e V>ulh,Vl} . 

Proof. See (2TJ Prop. 2.28, (7.29)]. □ 

From Lemma [38] and Proposition I35[ we can obtain additional estimates for the 
inverse Hodge operators; see [2TJ Prop. 2.29, (7.30)] for details. 

6.6. Proof Outline II: The Sharp Trace Estimate. The next batch of auxiliary 
estimates are the basic geometric Littlewood-Paley estimates and the corresponding 
Besov estimates on both N and the <S„'s. These were listed in [2TJ Sec. 7.3], while 
much of the background material is within [211 Ch. 2]. 

The main issue is that we are required to make use of Besov norms within the 
proof of Theoreml33"l Recall that in R 2 , although the Sobolev embedding H 1 ^ L°° 
fails, the corresponding Besov space B\ l5 on the other hand, does embed into L°°. 
In fact, we will need an analogous Besov embedding property on the <S„'s. This will 
be necessary in particular for the L^L^-type estimates for x and £. 

To begin with, one must make sense of Besov norms. This task was accomplished 
in [13] via the construction of a geometric tensorial Littlewood-Paley (abbreviated 
L-P) theory. In terms of our current setting, we can systematically define "L-P 
projections" Pfc, k € Z, which are used to dyadically decompose horizontal tensor 
fields on TV. These P/c's are defined explicitly using the heat flow operators intrinsic 
to the S v J s; see [13l Sec. 5] and [2.1, Sec. 2.2] for the exact definitions. 

One can proceed to show that these geometric L-P operators satisfy many of the 
same properties as the classical L-P projections on Euclidean spaces. 

Lemma 39. Let * e TTAf, k > 0, and < v < 5 . 



2 ^See also |21l Sec. 2.2] for additional remarks and clarifications on this construction. 
■^The idea of constructing Littlewood-Paley theories from the heat flow, or more generally 
from diffusion semigroups or martingales, originated from 1231 . 
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• The following estimate holds for each 1 < p < oo: 

(H2) \\Pknms v )Z\\m L z { s v )- 

• The following "finite band" estimates hold for each 1 < p < oo: 

(113) \\£p k n LHSv) = ii^*iilp (5 „) z 2 2fe 11*11^^) , 

• The following "finite band" estimates hold: 

(114) 1^*11^(5.5 + \\mn L *( Sv ) z 2" ii*iil 2(s „) , 

ll^llL 2(5 „)^2- fc ||y*|| L2(5u) . 

• If the (NO) Pt s ,A condition holds for sufficiently small 8q (with respect to 
Aq and the fundamental constants), then the following "weak Bernstein 
inequalities" hold for any 2 < p < oo and 1 < q < 2: 

(us) \\Pkn LHSv) z 2( i -t) fc \\nms v) > 11^*11^(5.) £ 2(t- i ) fc \\m L , [Sv) ■ 

Proof. See [13]. □ 

In particular, similar to the classical L-P projections, the Pfe's convert covariant 
derivatives into multiplication by a scaling factor (in the weaker sense of L p -norms). 
Several other basic properties for the Pfc's can be found in p~3l [21] . 

We can now define (geometric tensorial) Besov norms using the operators P k . 
The Bcsov-type norms we will need for Theorem 1331 are the following: 



k>0 

u*iiw) = E \\ p k*h*w + wnmM) ■ 

k>0 

For example, one can prove the following sharp Besov embedding estimate: 

(116) IM| £ « (A0 <II^IU (A0 + ||t;V|| Lri2(jV) , ^ec°°(AA). 

Remark. In fact, one of the main reasons why this line of reasoning is so highly 
technical is because the inequality (| 1 1 6[) and other sharp L°° -L 2 and L 2 -L l esti- 
mates cannot be established for nonscalar tensorial quantities. This is due to the 
lack of a satisfactory tensorial "Bochner estimate" controlling the second covariant 
derivative J7 2 in L 2 -norms by A. The reason for this is that we only have a very 
weak Lf- ' H^ 1 / 2 -type control on the Gauss curvatures of the S v 's. 

The primary tool we will require from our Besov norms is the following auxiliary 
estimate, which we refer to as the "sharp trace theorem" . 

Lemma 40. Assume the fN0) Pi j D) A condition holds for sufficiently small Sq and 
sufficiently large A (with respect to the fundamental constants) . If ty, P,Be TTAf, 
and if^J^ = Tt P + E, then the following "sharp trace" estimate holds: 

wn^L-w z ii*iiw(ao + ii p ii^ ( Ao + ii^iiw) . 

Proof. See [2TJ Lemma 7.22] and [HI [24] . □ 



fact, the first step is an application of 1 16t . 
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The goal will be to apply Lemma l40l to \ an d C The proof of Lemma l40l requires 
an extensive number of Littlewood-Paley decompositions, integrations by parts, and 
Besov estimates. This is a significant technical undertaking in its own right; in fact, 
all of [14] is dedicated to establishing these estimates. The main ingredient of this 
proof is the "bilinear trace theorems" , in which one bounds integrals along null 
generators of bilinear tensorial quantities; these were stated and proved for the 
vacuum case in [14] Thm. 4.9] and [24] Prop. 5.4]. For the E-S and E-M cases, the 
proofs of these estimates remain essentially unchanged; see [211 Lemma 7.20]. 

6.7. Proof Outline III: The Null Injectivity Radius. Another one of the 
most important auxiliary estimates is that of bounding the null injectivity radii 
from below. Such an estimate would state roughly the following: 

Lemma 41. Assume the bootstrap assumptions, i.e., that (N0) qj s ,A holds for any 
q G M + . If (5o is sufficiently small with respect to the fundamental constants, then 
the (INTJq^o.Ao condition also holds for every q € M + . In particular, this implies 
a uniform lower bound for the null injectivity radii i{q), q G M + . 

The proof of Lemma [41] (in the E-S and E-M cases) is essentially the same as 
the vacuum case, originally presented in [16 . Only a slight adaptation of this was 
required for our settings; this was outlined in [5TJ Sec. 7.4]. 

For completeness, we provide a short summary of the argument behind Lemma 
I4T1 We must control both the null conjugacy radius s(q) and the null cut locus 
radius [(g) of each point q G AI + . [3 First, we can control the null conjugacy radii 
by appealing to the following "breakdown principle" for regular null cones: 

A null conjugate point cannot occur along Af~ (p) if tr \ remains finite. 
This is true because uniform control on tr x implies control on the volume forms 
of the 5„'s, and hence the (w, cj)-parametrization of M~(p) can be shown to be 
nonsingular. As a result of this, we obtain the following: 

Lemma 42. If (N0) p ^ ,a holds, then s(p) > min([(p), Sq). 

Proof. See (2TJ Lemma 7.23]. □ 

It remains to control the null cut locus radii l(q), q G M + . The argument for 
this is identical to that found in [16] . The main ingredients are the following: 

• The regularity properties of the time foliation, established from the a priori 
estimates of Sectional In particular, this includes (|40|) and Proposition |2~T1 

• The following null cone comparison property: 

About any q G M , there exists a sufficiently large "almost Minkowski" 
coordinate system about q, within which the regular past null cone is 
contained within two Minkowski coordinate cones. 
For details, consult [2H Prop. 4.11, Prop. 4.12]. 
As a result of the above, we can make the following observations: 

• The iS„'s must be "comparable" to rescaled Euclidean spheres. 

• If the first point of intersection between two null generators occurs before 
null conjugate points, then the generators must intersect with angle 7r. For 
a proof of this property, see [TH] Lemma 3.1]. 



Recall that s(p) is the timespan before which there are no past null conjugate points of p, 
while [(p) is the timespan before which no two distinct null generators of Af~(p) can intersect. 



BREAKDOWN CRITERIA 



45 



• If two null generators of M~ (p) are opposite at p, then they cannot intersect 
until after a fixed timespan <5» , which depends on the regularity of the time 
foliation, and hence on the fundamental constants; see jTHl Prop. 3.4]. 

We can show that the above observations suffice to imply a lower bound on the 
[(g) 's. The idea is roughly the following. Using [16| Prop. 3.5], we can generate a 
pair of null geodesic segments that intersect at two points po, Qo such that: 

• The timespan between po and qo is at most that of any other pair of null 
geodesic segments that intersect twice. 

• The two segments connecting po and qo are opposite at both po and qo . 

Then, the preceding observations, with the help of Lemma 142) will suffice to control 
this minimal timespan between po and qo. As a result of this argument, then Lemma 
l4T1is properly established. For the details, see [161 Sec. 3]. 

As a result, we obtain a lower bound on the null injectivity radii as well as the 
(Nl) gi j 0l A conditions. This will be important in the remainder of the proof of 
Theorem 1331 since it is used to derive trace estimates, which bound quantities on 
the S v 's by quantities on the timeslices E T ; see [211 Prop. 7.3, Cor. 7.4, Cor. 7.5]. 

Remark. We remark that as stated in Lemma \41\ the lower bound for the null 
injectivity radius is contingent on the bootstrap assumptions. However, after com- 
pleting the bootstrap argument and hence proving that the (NO) q ,s ,A properties 
hold without preconditions, then we can apply Lemma [7JJ again. This controls the 
null injectivity radii without requiring the bootstrap assumptions a priori. 

6.8. Proof Outline IV: The Improved Estimates. With the necessary auxil- 
iary estimates established, we now turn to deriving improved estimates for the Ricci 
coefficients and the Gauss curvatures /C. This would complete our main bootstrap 
argument. The details sketched here are performed in detail in |21[ Ch. 8]. 

The Ricci coefficients are related to each other and to the spacetime curvature 
R via a family of evolutionary (along null generators) and elliptic identities, known 
as the structure equations. For the full list of structure equations, see [21[ Sec. 3.2]. 
These are adaptations of the analogous equations for vacuum spacetimes found in 
pj [TJ1 HH1 [2H ■ In particular, in the nonvacuum settings discussed here, there are 
numerous additional terms in the structure equations corresponding to the Ricci 
curvature; of course, these terms vanish in the vacuum setting. 

Our desired improved estimates will be a consequence of the forms of these 
structure equations. Let T denote one of the quantities bounded by the (N0)p !( 5 .A 
condition, i.e., T < Aq. Our goal, in general, will be to integrate the various 
structure equations in order to derive estimates roughly of the form 



where the constant of the inequality depends only on the fundamental constants. In 
general, the bound Sq^ 2 Aq is a consequence of terms in the structure equations that 
are quadratic (and occasionally cubic) in the Ricci coefficients, while "1" reflects 
those terms which can be controlled using a priori estimates. 

If we choose Ao to be sufficiently large with respect to the fundamental con- 
stants, and we then choose 5q to be sufficiently small with respect to Ao and the 
fundamental constants, then we obtain the schematic bound 



r < A 2 a 2 

1 ~ °0 ^0 



+ 1, 
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In particular, this means we have proved the (N0) Pj 5 0i a /2 condition. As a result, 
the bootstrap argument would be complete. We will now briefly describe how the 
above schematic bound can be obtained for the possible values of F. 
For example, consider the following structure equation: 

(117) f tp (tr x) = -\* (tr X ) 2 - \x? ~ VR*. 

This is the Raychaudhuri equation, which is perhaps the most well-known of the 
structure equations. By integrating a slight variant of the Raychaudhuri equation 
and taking an L^L^-norm of the result, we derive 



tf(trx)-- 



1. 



L~L? (AO 



The quantity A 2 , follows from the bootstrap assumption T < Ao, while the small 

1/2 

factor S appears because of the integrations involved. The Ricci curvature term 
in (|117p can be controlled using a priori estimates from the previous section. 

We can also control the quantity )Pt p [$(tox) — 2ip 1 ] in the L 2 (A/")-norm using 
(|117[) . Similar bounds for the L 2 -norms of ^ t p X an d ^t p C can be derived using the 
corresponding evolution equations for x an< i C 

At this point, we can also prove analogous estimates for the quantities in the left- 
hand side of (|110[) by integrating the appropriate evolutionary structure equations. 
These estimates will in fact function like the auxiliary estimates: 

• These bounds will be essential for completing the bootstrap argument. 

• Once the bootstrap argument is completed, and (|109l) is established, then 
(|110p will immediately follow from (|109[) and these estimates. 

Remark. The evolutionary structure equations are listed in [211 Prop. 3.8] . 

Next, we turn to the elliptic structure equations, cf. |21[ Prop. 3.9], from which 
we can control y/\ an d m the L 2 -norms. These equations are of the forms 

T>2X = R + f(trx) + B, 2?iC = R + /i + B, 

where R represents certain components of the spacetime Riemann and Ricci cur- 
vatures, and B represents lower-order terms. Since the right-hand sides of these 
elliptic equations can be easily controlled in the L 2 (A/")-norm, then direct applica- 
tions of Lemma 1551 yield the desired improved estimates for ^x an d ^C- These are 
once again of the familiar form V < Jq^A 2 , + 1. 
Furthermore, using the structure equation 



K = -- {tv X ) {^x) + \x ab X ab 



-R4343 + TT-Ra" 



we can proceed to derive the following improved estimates: 



-R 



K. 



1 



m 2 p 



1. 



1 



m 2 p 
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1. 



L 2 (A0 



The latter bound is straightforward, but the former requires much more effort due 
to the presence of the fractional Sobolev 7J _1 / 2 -norm. IfpJ Indeed, for the LfB^ 1 ! 2 - 
bound, we need a technical estimate for the commutator [A _1 / 2 ,y t ], where A -1 / 2 



33 See [13] and [21] Sec. 2.2] for the definition of the H s -norms. 
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is the fractional Sobolev operator on the iS w 's, defined in the standard fashion in 
terms of fi. For details, see [2TJ Lemma 8.8], as well as p~2l 124] . 

It remains only to establish improved estimates for the L^L^-norms of x an d 
£. As mentioned earlier, the strategy is to utilize Lemma 1401 In particular, we 
must show that fx and f( have decompositions of the form ftpP + E, where the 
quantities P and E have sufficient bounds, as specified in Lemma 1401 This is a 
priori not obvious; the argument in fact spans nearly all of [21[ Sec. 8.2-8.4]. 

The main ideas and steps for deriving these decompositions are the following: 

• First, since the structure equations only provide relations for T>2X arl d £>iCi 
not for "fx and J/^, we must apply the inverse Hodge operators, e.g., 

f x =fV^V 2 x = fV^R+yV^f (tr X ) +fV^B. 

• In order to generate the components ^t p P m the desired decomposition, 
we must take advantage of additional relations satisfied by R. For instance, 
for C, this is a "null Bianchi equation" of the form 

V\R=f tp R + D, 

where D denotes lower-order terms. From this, we see that we will also 
need to work with the "inverse" operator (T>\)~ 1 . 

• Next, combining the previous two steps, we see that we must commute 
y tp with f and the inverse Hodge operators. In particular, this requires 
multiple commutator estimates in both Besov and lower-order norms. 

• Unfortunately, there exist "bad" commutator terms which cannot be ad- 
equately controlled. In fact, these terms must themselves be further de- 
composed by adopting arguments analogous to the above, which yields an 
additional level of "good" and "bad" terms. This ultimately results in 
an infinite sequence of decompositions, which can be shown to converge. 
Only in this fashion can the "bad" commutator terms be treated. 

• After the infinite sequence of renormalizations, we finally obtain the desired 
decompositions for 'fx and fC,. We must also show that components of 
these decompositions can be controlled as required by Lemma l40l 

Once all of the above is complete, then Lemma [40] yields the desired improved 
estimates for both x and ( in the L^L^-norm. 

Remark. For the most part, the Ricci curvature terms in the structure equations 
(which are unique to the nonvacuum settings) are harmless, as they generally have 
better bounds than the Riemann curvature terms due to (|60p and Proposition \24\ 
The main exception to this is for the above decomposition argument, for which 
the higher- order Ricci curvature terms require their own unique decompositions. In 
particular, for the E-M case, we must make use of analogous null Bianchi equations 
for components of F. For details, see |211 Sec. 8.3]. 

As a result of all of the above, we obtain the improved estimates T < Ao/2, which 
proves the (N0) 9 ,s ,a /2 condition for all q G M + and hence completes the boot- 
strap process. By a standard continuity argument, we obtain that the (NO) ?j j ,A 
condition holds for all q G M + without the a priori bootstrap assumption. 

Finally, by Lemma [44] the conditions (Nl) g! 5 0i A also hold for all q G M + 
without prior assumptions. This implies the desired past null injectivity radius 



34 The "good" terms will be contribute to bothy tp P and E. 
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bounds, along with (jl09|) . The remaining estimates (|110|) now follow from ()109|) . 
as discussed before. This completes the proof of Theorem l33l 

7. Higher-Order Energy Estimates 

Now that we have established both a priori estimates and local null geometry es- 
timates, we can proceed to derive higher-order energy estimates and then complete 
the proof of Theorem [5] Our main goal is higher-order L 2 -estimates on k, 1Z, and 
the horizontal formulation of the matter field; these are the objects represented 
by the quantities &(t), $K(t), f(r) in pi|) - P6) . Like for the a priori estimates, 
we obtain these by first deriving higher-order energy inequalities for the spacetime 
quantities R and the matter field and then deriving elliptic estimates. 

Unlike for the a priori estimates, we will also require L°°-bounds for R and the 
matter field (D 2 (f> in the E-S case, DF in the E-M case) in order to complete the 
higher-order energy estimates. For this task, we will make use of the generalized 
representation formula for covariant tensorial wave equations discussed in [211 122] . 
This is an extension of the formula of S. Klainerman and I. Rodnianski in [15j . 
which was applied in an analogous fashion in the E-V case. However, the formula 
of [TS] fails in the E-M case due to additional first-order terms. 

7.1. Covariant Wave Relations. We have previously made use of the fact that 
the matter field (D(f> in the E-S case, and F in the E-M case) satisfies a covariant 
wave equation. We now list similar wave equations for R, DR, and higher deriva- 
tives of the matter field. We can then derive energy inequalities from these relations 
using the same EMT techniques as in Propositions [21] and [24l 

To begin with, by taking the divergence of the Bianchi identities (|85l) for R and 
then commuting derivatives, we obtain the covariant wave equation 

U g R a p lS = -D^Rp^s ~ Dp^R^S + Ra^ [Rlp^S + [ R ]^a-yS ■ 

Applying either (f25j) or (|26|). depending on the matter model at hand, and expand- 
ing the terms quadratic in i?, we obtain the nonlinear wave equation 

(118) D g R = d 2 Ric +R * R, 

(R * R) a p lS = 2la/3 {R a X Rxi3fs} + 2i? Q A p h l Rfj,XfS — 22l Q /3 {Ra^^RpiiSx} ■ 

The term d 2 Ric has the following expansions in the E-S and E-M cases, respectively: 
(dP Ric) = -22l Q/3 {D al <j>D 0S <j>} - 2l Q/3 2l 7 a {D^D^R^} , 

(d 2 Ric) = 2l Q £ {F^D^Fjs} + % s {F^D s „F a0 } - Vl a p% s {D a F^D 5 F^} 
+ D^F a0 D^F lS + 2t Q ^2t 75 {F x ^F^R aXps + 
- ^ a /3%s {g ai (F^DpsF^x + DpF^DsF^)} , 

Recall that 21, defined in (J7J, denotes anti-symmetrizations of the indices. 

In the E-S case, by differentiating (|86|) and commuting derivatives, we can derive 
a wave equation for D 2 (p. More specifically, we can compute 

n g D a0 <l> = £>V^ + (Rp\ a D x ^) 

= D a U g D 4> + R^ a D Xfj( f> + R x " a D„ x( f> 

+ D^R^ aP x D x <j) + R^D^cj). 
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By ([25]) and (J86j>, then 

(119) = 2D X (f>D x( (>D a f j( f> + D a <j)D X (j>D P x<i> 

+ D P 4>D x 4>D aX <t> - 2R a »p X D flX <f>. 
By a similar process, in the E-M case, we have the wave relation 

(120) D g D^F a p = 2F x »D 7 R aX ^ + 2lo^ {Ra^^D^Fpx} + 2R a Xfi pDjFxfj. 

+ 2F a »Fp X D 1 F tlX + 22l Q(9 {F^ x F a ^F 0X } 

- F^FcpD^x ~ ^F» x F^D 7 F a p 

- F x ^F ltl D x F aP - 2t Q/3 {F a x F^D^F px } 

- 2l Q/3 F a A 2t^ A (V^£ A F 7M - ^gp^F^ D X F^ 

In addition to the above, we will need covariant wave equations for one higher 
derivative of R, (f>, and F. Fortunately, we will require only the schematic forms, 
not the exact equations. In the E-S case, we have 

(121) U g DR ^ D 2 (j) ■ D 3 cj) + D(j> ■ D 2 cj) ■ R + (Dcj)) 2 DR + R- DR, 

(122) D g D 3 <j) = (Defy 2 ■ D z (f) + D(j> ■ (D 2 (f>) 2 + R ■ D 3 (f) + DR ■ D 2 (f>. 
Similarly, in the E-M case, 

(123) U g DR ^ F ■ D 3 F + DF ■ D 2 F + F ■ DF ■ R + F 2 ■ DR + R ■ DR, 

(124) D g D 2 F F ■ D 2 R + DF ■ DR + R- D 2 F + F ■ (DF) 2 + F 2 ■ D 2 F. 

7.2. Energy Inequalities. With the wave relations (|118[) - (|124j) in hand, we can 
now apply general EMT methods to derive preliminary energy inequalities for the 
curvature and the matter field. Like for the a priori estimates of Proposition [2~T1 
we must bound the curvature and the matter field concurrently. 
For convenience, in the E-S case, we define the quantities 



£ 2 (r) = [I-DJZII^) + ||£> 3 C=(£ T ) . S * W = \\ D2R \\ 2 m^) + \\ D4 4U^) 



while in the E-M case, we define the analogous quantities 



£ 2 (r) = WDRf^ + \\D 2 Ff L2{ ^ } , £» (r) = ^R^^ + U^ffc^ 



Proposition 43. For any r < T\ < r < t\, the energy inequalities 

£ 2 (r) < S 2 [n) + f (l + WRWU^ + ll^ll^^,)) dr', 

J Tl 

£ 2 (r) < £ 2 (n) + f (l+ \\R\\l°o { -z Tl ) + \\DF\\l~ ( z T ,)) dr' 
hold in the E-S and E-M cases, respectively. 

Proof. In the E-S case, applying (|83p with $ = R along with (|118l) yields 

\\DR\\l H x T) £ \\DR\\l^ T1 ) + r f (\DR\ 2 + \R\' + \D 2 <J>\ 4 + \R\ 2 ) dr. 
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Similarly, applying ((83]) with $ = D 2 (f> along with (|119[) yields 

w^fm^ s ii^ 3 ^iIl 2(Sti ) + / r ^ (l^>| 2 + l^l 2 l^>| 2 + l^>| 2 ) 

In these inequalities, we have used the a priori uniform bounds on k, n, and Dip. 
Summing the above and applying Proposition 1211 yields 

£ 2 (r) < £ 2 (n) + (l + ||i?||^ (ST , ) + \\D^t L ^ ri) ) dr' + £ £ 2 (r') dr' . 

Applying Gronwall's inequality results in the desired bound. 

The E-M case proceeds analogously; apply ([83]) with $ = R and $ = DF. □ 

The process behind the proof of Proposition [43] can be repeated for the wave 
equations (|121j) - (|124[) in order to obtain the following: 

Proposition 44. For any tq < t± < r < t\, the energy inequalities 

£ 3 (r) < £ 3 (n) + f (l + 11^11^-^) + \\D 2 4 2 L ) [1 + £ 2 (r')] dr' , 

& (r) < £ 3 (n) + [(} + Pll W T ,) + \\DF\\ 2 l^ t ,)) [l + ? V)] dr' 
hold in the E-S and E-M cases, respectively. 

Proof. In the E-S case, apply ([83]), (fT2l]> . and (jT22]) with $ = DR and $ = L> 3 0. 
In the E-M case, apply ((HH) , ([123]) . and (|T24| with $ = Di? and $ = D 2 ^. □ 

Remark. FFe can obtain identical energy bounds for DR and D 2 R, as well as DF, 
D 2 F, and D 3 F in the E-M case, using the generalized Maxwell EMTs and (|84[). 

Lastly, we will need a local variant of Proposition [43] involving flux bounds on 
regular past null cones. This can be derived in the same manner that Proposition 
[24] was obtained from Proposition [21] Fix p £ M + and a time value tq < t < t(j>), 
with t{p) — r < lip); we normalize and foliate N~(p) using T\ p and t p . 

Define the sets Af T , I T , and as in Section [5751 and define the flux densities 

p 2 ip, t) = -Q w [R] (T, L) - Q w [D 2 4>] (T, L) , 
P 2 (p, r) = -Q w [R] (T, L) - Q w [DF] (T, L) , 
in the E-S and E-M settings, respectively. We also define the corresponding flux 

J 72 ip; t) = / p 2 ip, t) . 

Recall that the notation Qu,[J7] was defined in Section [5731 

Proposition 45. The following hold in the E-S and E-M settings, respectively: 
F 2 ip;r)<£ 2 iT) + J (1 + PHioc^,) + ||£ 2 Coo (IV) ) dr', 
(p; r) < £ 2 (r) + J (l + p]]^,) + 11^11^^,)) dr' . 
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Proof. For the E-S case, we apply ([88} with $ = R and $ = D 2 (j) to obtain 

J (\DR\ 2 + \Rf +\D 2 ^ + \R\ 2 )dT', 
-J Q w [D 2 cp](T,L)<\\D^f L ^ T) 

(\D 3 <f>\ 2 + \R\ 2 \D 2 ^\ 2 + \D 2 ^ dr'. 



Summing the above bounds and applying Proposition [211 we obtain 



f 2 (p;t) <£ 2 (r) + 



Hp) r 



l + ll^l 



i°°(S T /) 



dr'. 



The proof is completed by applying Proposition H51 

The E-M case is analogous, except we apply (|88p with $ = i? and $ = DF. □ 

7.3. The Generalized Kirchhoff-Sobolev Parametrix. From Propositions [43l 
and 2H we see that in order to control the higher-order energy quantities £ 2 (r) 
and £ 3 {t), uniformly for to < r < t\, we must first obtain a bound for the L°°- 
norms of R and the matter field: D 2 (j) in the E-S case, and DF in the E-M case. 
As mentioned before, we will accomplish this by using the wave equations (|118D . 
(|119[) . and (|120[) . along with the generalized Kirchhoff-Sobolev parametrix of [22] . 

In the vacuum analogue of [17], one needed an L°°-bound for R, which was 
obtained using the Kirchhoff-Sobolev parametrix of |15j . Although this parametrix 
also suffices in the E-S setting, it fails in the E-M case due to first-order terms 
present in (|118[) and (|120|) . Such terms must be handled differently than allowed 
in [15] ; this is the primary motivation behind the generalized formula of |22j . 

Consider the scalar wave equation in the Minkowski spacetime R 1+3 : 

(125) n<f> = ip, 0| t=o ="o, d ff>\ t =o = a i- 

From standard theory, cf. [TQl Sec. 2.4], at a point (t,x) S (0, oo) x R 3 , we can 
express (f>(t,x) explicitly in terms of ip, cto, and a\. Moreover, the above can be 
written as an integral along the past null cone segment N from the vertex (i, x) to 
the timeslice t = 0, along with an "initial data" integral on the set N n {t = 0}. 

The goals of both [21] and its predecessor [T5] are to provide a local first-order 
extension of the standard formula to arbitrary curved spacetimes. In particular, 
the parametrices of [15, 22 enjoy the following features: 

• These formulas treat tensorial wave equations in a completely covariant 
fashion, without first expressing them as scalar equations. 

• The parametrices are supported entirely on past null cones. 

• The parametrices are valid only on regular (i.e., smooth) portions of past 
null cones. In the setting of this paper, this means that the formulas are 
only applicable up to the null injectivity radius. 

• The parametrices contain "error terms" , expressed as integrals along regular 
past null cones. These are a result of the nontrivial geometry. 

• The parametrices depend only on quantities defined on the past null cones, 
i.e., it is independent of the extensions of such quantities off the cones. 
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The parametrices can be systematically generalized to covariant wave equa- 
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tions on arbitrary vector bundles over an arbitrary curved spacctime. 
The parametrix of 22 differs from that of [15] in the following respects: 

• The formula of |22j directly handles systems of tensor wave equations with 
additional first-order terms. The first-order terms are treated by altering 
the transport equation associated with the parametrix. 

• While the proofs in |15j make heavy use of distributions on manifolds in an 
informal fashion, the proofs in |22j instead remain entirely at the level of 
rigorous calculus operations on null cones. 

• The computations in [T5] generated terms not supported on the null cone, 
which must all be meticulously cancelled in later steps. In [22], on the other 
hand, the entire derivation was performed on the null cone. More specifi- 
cally, the extra terms in [15] could be avoided altogether if one integrates 
by parts only the derivatives tangential to the null cone. 

• As a consequence of the above modifications, one can slightly weaken the 
assumptions listed in |15] for the parametrix to be valid. I " 

For more details on the preceding discussions, see [22], or [21] Ch. 5]. 

We now state a special case of the main result of [22], adapted to the setting 
of this paper. More specifically, we state the representation formula only in the 
case of time foliations of the null cone, while [22] also considered arbitrary foliating 
functions. For the general version, see [32] Thm. 7]. 

Theorem 46. Assume the following: 

• Let n be a positive integer, let r (1) , . . . , r ( ™' be nonnegative integers, and 
suppose for each 1 < m, c < n, we have defined tensor fields 

$(m) £ YT r(m) M, * (m) 6 TT rim) M, P (mc) € rT 1+r<m)+r<c) M. 

• Suppose the <J> (m) 's, <3> (m) 's, and P< mc > 's satisfy the system 

n 

(126) □ s $ (m) / + ^P <mc V/ J D Al $ (c) ,/ = * (m) /, l<m<n 

c=l 

of tensor wave equations, where I and J in (|126[) are collections of r <m) and 
r <c) spacetime indices, respectively. 

• Fix p 6 M, and suppose Af~ (p) is normalized and foliated by T\ p and t p . 

• Let vq be a constant such that < vq < i(p). 

• For each 1 < m < n, we define the extrinsic tensor fields 



B im) e rT - <m V- (p) , A lm) = t~ G TT r{m) N~ (p) , 



along with a tensor J <m) of rank r <m) at p, such that the system 

i) 
2 



(127) f tp B^ = -\ 



p 



B (m)/ + |^P ( ™'4/S <c)J ; l<m<n 

c=l 



of transport equations is satisfied, along with the initial conditions 
(128) B (m} \ p = J (m) ; l<m<n. 



^In this extended setting, one also needs a bundle metric and a compatible connection. 
^In fact, this weakening of the assumptions applies to both the formulas of 1151 and 1221 . since 
the results in 1 151 are strictly special cases of those in 1221 . 
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• Define the "error coefficients" 

„<«») G T fr^+r^ N - ^ } \< m ,c<n 

by the formulas 

(129) v^/ = -TP^a/ + IfiP^SJ 1 + CP^a./ + \ (trx) P^Uj 1 



n 

+ i (tr X ) P (em) 3 / + \ E P^^P^hK 1 - 

d=l 

Then, we have the representation formula 

n 

(130) 4tt • n (p) • E Jlm)I & m) i\ P = $ (p; «b) + e 1 (p; «d) + (p; «d) + 3 (p; «o) , 

where: 

• The '[fundamental solution term' 7 $(p;vo) is given by 

n „ 

(131) d(p;v ) = -J2 / A^^^j. 

m=1 JM-(p;v ) 

• The "principal error terms" € 1 (p;v n ) are given by 

n 

(132) £ 1 (p;vo) = -Y, / TA {m)I fa^ m) i 

TO =i (p;«o) 

n 

+ E / (C a -^ a )f a ^ (m)7 $ (m) J- 



i N (p;vo) 

The remaining "error terms" <£ 2 (p;vo) are given by 



n ,. 

(133) £ 2 (p;v ) = E / ^ ' A {m)I <& {m) i 

m=1 JAr-(p;v ) 

+ ^E / A ™ lR " [ $(m) ]/ 

n „ 

- E / p"" 1 ^""*'"', 

n „ 

+ E / V Cem) J 7. i4 (e)J $ (m) J _ 
m,c=l"' A '' _ (P;' u ") 

• TTie "initial value terms" 3(p;v ) are given by 

(134) 3(p;v ) = -^f2 [ (trx)^ ,7 * (m) /-E / 4 (m)I I>3* 

m — 1 JS V0 m=1 J5„ 
n . 

-5 E / p (era> 3j^ ce)J * (m) /. 

m,c=l "^"o 

-ffere, we /iawe indexed with respect to arbitrary null frames L,L, ei,e2 adapted to 
the t p -foliation. The capital letters I, J refer to collections of extrinsic indices. The 
symbols \, Xj C; aKC ^ »? refer to i/ie fficci coefficients of N~{p). 
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Remark. It is easy to see that Theorem \46\ is a special case of [22\ Thm. 7] . 

The function t p is clearly a foliating function in the sense of [221 Sec. 2.2], and 
the associated null lapse $ has initial value n(p) at p. Moreover, both Af~(p) and 
J\f~(p]Vo), as given in Theorem \46\ are regular portions of the full null cone N~(p). 

Remark. Although the representation formula was stated in (|130|) - (|134j) in index 
notation, this was done only as a matter of convenience. It is easy to see that these 
expressions can in fact be described invariantly. 

Remark. In particular, we can use (|130j) to examine the value of any $ (m) |p indi- 
vidually by setting J <c) = for all c ^ m. 

Lastly, we note that the ^4 (m) 's in Theorem Flol satisfy the transport equations 

i i " 

(135) f L A^ I = --(t rx )A^ I + -Y / P <cm) 4J I A Mj , l<m<n. 

c=l 

7.4. Applying the Parametrix. We now describe how Theorem P||>] is applied. 
Let < 5q < 1 be sufficiently small such that Theorem 1551 is satisfied. |f] Fix a point 
p € M_|_, fix another constant < S < min((5o, t(p) — To), and define N = Af~(p; S), 
normalized and foliated by T\ p and t p , as usual. 

In the E-S case, we wish to find L°°-bounds for both R and D 2 (j) at p. As a 
result, appropriating the notations of Theorem l46l we set n = 2, = R, and 
$ (2) = D 2 <j>. For the corresponding system of wave equations comprising (|126[) . we 
take (fTTH)) and (Ti~T9|) . Moreover, the right-hand sides of pTg)) and (fTT9|) determine 
the v& (m) 's and P <cm) 's. In particular, in the E-S case, the wave equations contain 
no first-order terms, that is, the P< cm >'s vanish entirely. [3 

The E-M case is analogous, except that we must obtain L°°-bounds for R and 
DF. We apply Theorem by setting n = 2, = R, and $ (2) = DF, and we 
adopt (|118|) and (|120|) as the system of wave equations. In contrast to the E-S 
setting, the first-order coefficients p( cm > are no longer trivial. In particular, we see 
from the right-hand sides of (fTTH]) and (fT20|) that P (11) and P (22) vanish, while P (12) 
and P (21) are sums of terms, each of which can be expressed as tensor products and 
contractions of F with instances of g. More specifically, the P <12) and P (21) terms 
arise from the first, second, and seventh terms of (|118[) in the E-M expansion of 
d 2 Ric, and from the first term of (|120[) . The remaining terms on the right-hand 
sides of (|118l) and (I120p comprise the ^"''s. 

Now that we have determined the P (cm) 's, we can determine their a priori bounds 
on J\f. By the breakdown criterion (|40j) , we have the uniform bounds 

(136) \\P (cm) \\ L ~ W Zh l<m,c<2 

in both the E-S and E-M settings. Furthermore, it follows from Propositions |2~41 
and [25] that in both the E-S and E-M cases, 

(137) |F^ m) |L 2(A0 + IIF^-il^AO ~ 1; 1 < m ,c<2. 
Of course, in the E-S case, the bounds (|136[) and (|137|) hold trivially. 



■^More explicitly, we assume that given any q g M+, then i(q) > <5i = min(<5o, t(g) — to): an d 
the estimates H109I I and (1 1 101 ) hold on Af q = M~{q\ <5i). 

"^Consequently, the parametrix of 1151 suffices in the E-S setting. 
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Finally, we choose arbitrary tensors J (1) and J (2) at p, of the same ranks as 
and <J> <2) , respectively. For convenience, we adopt the abbreviation 

|J| = |JW| + |JW|. 

Applying Theorem HHl we obtain the expansion (|130[) for the quantity 

47r-n(p)-(jW z ^ 1 > J | p + J< 9 > J $Wj|p). 

Each term of this expansion is either an integral over J\f (i.e., the terms $(p;8), 
(S (p;5), and <£ 2 (p;8)) or an "initial value" integral over the t p -level set Ss (i.e., 
the terms 3(p; 5)). Thus, in order to bound and $ (2) at p, as desired, we must 
bound each of the integral terms mentioned above. 

Before we can accomplish this, though, we must first control the fields B (m) and 
A (m> associated with the system of transport equations (|127|) . (|128[) . and (|135|) . This 
is the content of the subsequent proposition, which is the E-S and E-M analogue of 
the estimates [111 Prop. 7.1, Prop. 7.3] in the vacuum case. 

Proposition 47. In both the E-S and E-M settings, we have 

2 2 
m—1 m—1 

as long as So is sufficiently small with respect to the fundamental constants. 
Proof. By standard calculus computations analogous to those of Proposition [H 
ft P \B^\ 2 < \B™\ \f tp B^\ + \f tp h\ \B^\ 2 



< I5 (m, l 



\I\\B (m> \+J2\P^\\B^ 



\B {m) \ 2 . 



where I = d(tr x) ~ %t p 1 and h is the induced Riemannian metric on M. Here, we 
have applied (I3T))) and (I127p . Summing over m and applying (jllOl) and (I136[) . then 

jtJi2\ Blm) \ 2 )z( 1 + \ I \ + E i p(cm) i) Eis (m) i 2 <E^ <m, i 2 - 

\m— 1 / \ c,m— 1 / m—1 m—1 

Applying Gronwall's inequality to the above while taking into account the initial 
conditions (|128j) . we obtain the first desired estimate. 
Next, we define 

U im) = t 2 p -fA {m) e TT}T Am) N- (p) , |V| = max (\f A (1> | , \f A {2) \) . 

Differcntating the transport equations (|135[) and commuting derivatives (see the 
commutation formula |21[ Prop. 3.12]), we obtain the transport equations 

2 



C K U^ aC 



(138) f tp U {m) aK = -I ■ U {m) aK - tiXaU^hx + ^E^^ 4 

1 c=l 

- h p f a I ■ B {m) k + t v ti Y R^ ia B^ K , 
z ^ 



c=l 
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where K — (71, . . . , 7 r (m) ) and C denote collections of r <m) and r (c) extrinsic indices, 
respectively, and where denotes K except with the i-th index replaced by /1. In 
addition, we index both the horizontal and the extrinsic components of the U im) 's 
using local Fermi transported null frames. [3 If we also (locally) define 

C = \D L ei\ + \D L e 2 \ + \D L L\ + \D L L\ , 

then C < 1 77 1 < 1 by Proposition [3JJ and 

(139) \ft p (U^ aK )\ < \J tp U^\+C\U^\ < \f tpU ^\ + \U^\, 

where the left-hand side denotes applied to the scalar quantity U {m) a K- 
With consideration of the coefficients present in (|138[) . we define 

2 

Qi = |i| + lxl+ E l p(cm, l' 

c,m— 1 

2 2 

Q 2 = \fi\+m\ 2 ip ( ™ ) i-Q™[i?](T,L)+ e i?p 



(cm) I 

c,m— 1 c,m— 1 



In addition, by Proposition EMI (I109p . Proposition [371 (I136p . and (|137p . we obtain 

(140) HQlll^(AT)^l. IIQ2|| i2(A0 < 1. 

We now integrate (|138[) along each null generator of M. By also noting Corollary 
[26] and ([139]) . then we obtain the inequality 



v 2 \V\\ M < fw 2 (1 + Qi) 

« 



+ l J l / w fia| (lBlW )dti> 

« 



< 



w 4 \v\' 







+ |J| / w QaL.^dw, 



(I 



for each < u < 5 and u; G § 2 , where we also applied the preceding bound 
\B {m) \ < \ J\ and the estimate (|140p for Qi. Dividing both sides of the above by u, 
taking L^L^-norms of the resulting inequality (cf. (jTUTJ) ) , and then recalling the 
norm comparisons of Proposition (M[ we obtain 



5 

dw } dutdv 



\\\V\\\l^)< I I ( / v?\V\ 



JS 2 vo 



(to, a;) 



|J ' 2 LJo V Jo wQa ^ du ' ] dv ' U ' 

< JT HIT -II - I T|2 / / „,2 



< 



W\\\ 2 L , W + \J 



2 



where we applied Proposition [M] to the y-term, and where we applied Hardy's 
inequality and (|140p to the Q2-term. Shrinking 5q if necessary (still depending only 
on the fundamental constants), then || |V| ||i»(j\n < |J|, completing the proof. □ 



:S!) 



By this, we mean null frames L, L, ei, e2 satisfying the conditions ^i,e a = 0. 



BREAKDOWN CRITERIA 57 



7.5. The Uniform Bounds. We described in Section 17.41 how Theorem [46] is 
applied to the breakdown problem. With the preliminary work in place, we must 
now bound each of the resulting terms $(p,S), € x (p, 8), <B 2 (p,S), and 3(p,5). For 
convenience, let 8 denote D 2 <f) in the E-S case and DF in the E-M case. In addition, 
we define for every tq <r<t\ the uniform norm 

S{r) = \\R\\ L ^ T) + \\Q\\ L ^y 

We start with the principal "fundamental solution" term $(p; S), given by (|131[) . 
First of all, we decompose the integrands A {m)I ^> {m> j using local null frames. Since 
tpA (m> — B (m) is uniformly bounded by Proposition!!?! then the main challenge is to 
bound the integral of every component t" 1 ^/'" 1 '/ along TV. The primary observation 
is the following, which is also the cornerstone of [HI HZ]- From (jll8p - (|120[) , we see 
that each term of any component $ <m) / is of the schematic form W 2 ■ Q or 
where ||W||£°°(.A0 ^ 1> an d where Q denotes a tensor field of the form R or 0. 

The first form W 2 ■ Q can be sufficiently controlled as follows: 

/ \A^\\W\ 2 \Q\<\J\ [ t; 1 |Q| < a 2 |J| sup S(t). 

JM JM t(p)S<T<t(p) 

Note in particular that we have applied Proposition 27] to control the A (m), s. The 
terms of type Q 2 require a more detailed null frame decomposition. The main 
observation is the following: for each component of every Q 2 -term, one of the Q- 
factors is bounded by the flux density [p 1 (p,t(p) — which can be controlled 

using Proposition [Ml We shall further elaborate on this point later in this section. 
Assuming for now the above observation, then 

i 

/ A^-Q 2 < \J\[T 1 {p-t{p)-8)Y ( f t- 2 \Q\ 2 \ <6*\J\ sup S (r) 
JM \JM J t{p)-S<T<t( P ) 

by Proposition [24] Combining the above, we obtain the estimate 
(141) 3(p;S)<6?\J\ sup S(t)<8$\J\ sup 5(t). 

t{p)-8<T<t(p) t(p)-S<T<t(p) 

Remark. Recall that So > is the constant, depending only on the fundamental 
constants, for which Theorem \3S\ holds, while < i5 < mm(<5o, t(p) — To). 

We now discuss handling the Q 2 -terms in further detail. The process is analogous 
to the E-V case of |17) , except we must consider several additional terms involving 
the matter field. Recall that the exact null frame components controlled by the 
flux density p 1 {p,t(p) — S) were determined in detail in Proposition 1251 and in the 
subsequent discussions. More specifically, the only components of R and that 
are not controlled are those of the forms i?3 a 3b, D 33 (f> (in the E-S case), and D 3 F 3a 
(in the E-M case). As a result, we must determine that for every <5 2 -term in the 
^ <m) 's, at least one of the Q's is not an "invalid" component. 

We give a few examples here demonstrating the reasoning described above: 
• Consider the term / Q ^ 7< 5 = % a p[D ai ({)Dp$(j)\ in the expansion of (]118[) in the 
E-S case. By the above characterizations, the only invalid component, for 
which neither factor is controlled, is D 33 (j)D 33 (j). To obtain this component, 

"Recall that in the E-M case, the first-order terms (corresponding to the p( cm ) 's) are omitted 
from the 's. In fact, the reason that the original parametrix of |15l fails here is precisely 

because these first-order terms are not of the above forms W 2 ■ Q or Q 2 . 
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we require a = /3 = 7 = <5 = 3; but, by the antisymmetry between a and /?, 
the component 13333 vanishes. Consequently, any component of I a p-yS can 
be bounded by -Q w [D<f>](T,L) ■ \D 2 cj)\ < \p r \p,t{p) - 5)] 1 ' 2 ■ \D 2 cj)\. 

• The term I a p-yS — D^FapD^F-yg in the E-M expansion of (|118[) can be 
handled similarly. Because of the contraction involving fi and the properties 
of null frames, for each term of the summation, one of the /i's must not be 
"3". As a result, one DF factor can always be bounded by —Q W [F](T,L). 

• The terms R a x p^R^x-yS and ^a^Ru^^RptiSx] in ()118|) correspond with 
the i? 2 -terms present in the E-V setting in |17| . They can be handled using 
signature considerations, as in |171 Sec. 5.4], or directly using the same 
reasoning as for the previous two examples. More explicitly, we see that at 
least one of the R's is controlled by —Q m [R](T, L). 

• For the term R a ^ p x D in (|119p . due to the contractions, for each term 
in this summation, the "D 2 4>" factor can be bounded by —Q w [D(j)](T, L) 
unless /i = A = 3 in D 2 <f>. If this is true, however, then fi = A = 4 in the 
"i?" factor, which can hence be bounded by — Q m [R](T, L). 

The remaining Q 2 -terms in (|118p - (|120[) can be handled using similar reasoning. 

The next task is to bound the "error terms" ^{p^S) and <E 2 (p;5), defined in 
(fH21) and (fH3"]) . First, by Propositions 05] and S3 we have 



\fA (m) I Iy7$(" 



L 2 (J\f) 



< 



£ 2 {t{p)~5) 2 +62 sup S(t) 

t(p)-S<T<t(p) 



Similarly, by (|109|) and Proposition 1471 
|C-2| |y^ <m) | |$ (m) 

<s\J\ 

Therefore, we obtain the bound 



sup 

t(p)-s<r<t(p) 



\\L 2 (M) 

S(r). 



(142) 



< 



[£ 2 (t(p)-6)}- 2 +Sg\J\ sup S(t). 

t(p)-S<T<t(p) 



Next, for the first term of (I133[) . 

N|A«||$«|<| 



< 8 2 \J\ sup 

t( P )-S<T<t(p) 



L 2 (J\T) 
S{r) 



where we appealed to (|110|) to bound /i. For the second term, by expanding 
-R43[$ (m) ] using null frames, we see that each term of this expansion is bounded 
by -Q m [R](T,L) ■ |$ (m »|. As a result, 



A'' 



|^ m '||i? 43 [$<"»] I <^ 



J| sup S(t). 

t(p)-6<r<t(p) 



By (|136[) and Proposition Wf\ we can bound the third term: 



7A M < ||P ( ' 



in 



1 1/ 2 (AT) 



\l 2 (N) 
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<<J*|J| sup 5(r). 

t(p)S<T<t(p) 

For the last term, we will need the following simple estimate: 

Lemma 48. The coefficients v (cm) , as defined in (j!29[) . satisfy 

W^WlhM)^ 1 ' c, me {1,2}. 

Proof. We bound each term on the right-hand side of (|129[) . The first two terms 
are trivially bounded using (lf37j) ; by Theorem [33l Proposition [371 and (|136[) . the 

1/2 

third, fourth, and fifth terms are bounded by S Q . The final term on the right-hand 
side of (|129[) is trivially bounded using (|136l) . □ 

Finally, applying Lemma l48l then 

1^1 \ A ^\ < \\a^\\ L2W \\*™\\ L -w 



Af 



<d*|J| sup 5(r) 

t(p)-8<T<t(p) 



As a result, 



(143) <E 2 (p;6)<5*\J\ sup <S(r). 

*(p)->5<T<t(p) 

It remains only to bound the initial value terms 3(p; 5) in (|134[) . The first term 
on the right-hand side of (jl34[) can be handled using Proposition [34] and [37] 



f {tvx)\A^\\^\<S- 2 \J\ ( |$ M |<|J|||$ M || L «, (St(p) _ 
Applying (|T2"]) and Proposition |2~T1 then the above is bounded by 

\J\ [l + £ 2 (t(p)-S)+£ 3 (t(p) -*)]*. 
The third term of (|134|) can be bounded similarly: 

|pc«0| [AW I |$ <m) | < <5 1 J\ [1 + £ 2 (t 0) - 5) + £ 3 (t (p) - 5)] 



'5 5 

Lastly, for the remaining term of (|134[) . we will need the following trace estimate: 
Lemma 49. If Z E TTN~(p) and W S TTM are of the same rank, then 



Z L Wi 



ll^ll^( 55 )(l|v^L 2(St(p) _ 6) + II^IL 2(StCp) _ 5 )) 



Proof. See [21] Prop. 7.3]. The proof involves constructing a radial foliation of 
^t(p)-i5 using the level sets of a one-parameter family of null cones. This relies 
heavily on Theorem 1331 in particular the null injectivity radius bounds. □ 

We now apply Lemma [49] with Z = L ® A (m) and W — D<f> (m) . Since Z, as 
defined above, satisfies ||Z|| L oo( 5a ) < <S _1 | J|, then 



J L a A {m)I D a ¥ m) ! <6? \J\ (\\D 2 & m) \\ L2 ^ t(p} 5) + P$ (m) | 

<^|J| [£» {t{p) _ S) + g3 {t{p) _ S )]K 



L2 ( E t(p)-s) 
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As a result, we obtain 

(144) 3 (p; 8) < \J\ [1 + £ 2 (t (p) - 8) + £ 3 (t (p) -*)]*. 
Combining ([130} and (fT4T1) - ([T44|) . we obtain 

<4|J| sup S(t) 

t(p)-5<T<t{p) 

+ \J\ [l + £ 2 (t(p)-8) + £ 3 (t(p)-S)]K 
If we consider all possible values of J (1) and J (2) , then we have 

(145) \R\\ p +\e\\ p <sl sup S(r)+[l + £ 2 (t(p)-S)+£ 3 (t(p)-S)]K 

t{p)- s < T <t(p) 

Finally, if we fix r such that tq < t < t\ and r — 8 > To, and we apply (|145j) to 
every p £ M + satisfying r — 8 < t(p) < r, then we obtain 

sup S(t')<sI sup S(t')+ [1 + £ 2 (t-8) + £ 3 {t-8)]K 

t — 5<t'<t t—5<t'<t 

By taking Sq sufficiently small, we have proven the following: 

Proposition 50. Fix to < t < t\ . Suppose So > is sufficiently small, depending 
only on the fundamental constants, and let < 6 < So such that t — 5 > To. Then, 

sup S (r') 2 < 1 + £ 2 (t - S) + £ 3 (r - S) . 

t — S<t'<t 

7.6. Completion of the Proof. With the uniform bounds of Proposition [5U1 in 
place, we can now embark on the final stretch of the proof of the main theorem. We 
begin by applying Proposition [SO] to the higher-order energy estimates of Proposi- 
tions [43] and |44j This yields higher-order energy estimates. 

Proposition 51. For every To < t < t\, we have 

£ 2 (t)+£ 3 (t) <1. 

Proof. Let To < r < t\ and < S < So such that r — 5 > to, where So was 
determined in Proposition [3DJ Applying Proposition [SD] to Proposition 1431 yields 

£ 2 (t) <£ 2 (t-5) + 5[1+£ 2 (t-S)+ £ 3 (r - S)] 
<1 + £ 2 (t-S) + £ 3 (t-S). 
Combining the above with Propositions and [5D1 then 

£ 3 (r) < £ 3 (r - 6) + 5 [l + £ 2 (r - 5) + £ 3 (r - S)] 2 
<1+ [£ 2 {t-5) + £ 3 {t-S)] 2 . 
As a result, we have derived 

(146) £ 2 (r) + £ 3 (r) < 1 + [£ 2 (r - 5) + £ 3 (r - 5)] 2 . 

We can now iterate (|146[) . First, we let 8' = min(<5o, t\ — to), and we apply (|146p 
with < S < S' and t < tq + 8' to obtain 

sup [£ 2 {t') + £ 3 (t')] < 1. 

to<t'<to+i5' 



5^ J^'^'jI 

l<a<2 
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If 6' = t\ — to (i.e., t\ — to < So), then the proof is complete. On the other hand, if 
5' = So < t\ — to, then (|146[) also yields 

£ 2 (t + 6 q )+£ 3 (t q + 5 ) < 1. 

We can now repeat the above procedure, but with To replaced by To + <5o- By this 
process, we can reach the breakdown time t\ within a finite number of iterations. 
This implies the desired result and completes the proof. □ 

Now that the higher-order spacetime quantities are controlled, we must do the 
same for the corresponding horizontal quantities. Similar to the a priori estimates, 
we accomplish this by relating the corresponding spacetime and horizontal quanti- 
ties to each other and then applying elliptic estimates. The process is analogous to 
that of Sections 15.61 and [5 . 7\ except the computations are somewhat more involved. 
Consequently, we only list the final results and omit the proofs. 

Proposition 52. For any To <r<t\, the following inequalities hold: 

$K(T)+£(T)+f(Y) <1. 

Proof. This is a consequence of [STJ Lemma 9.4] and Proposition 15 II □ 

Remark. The quantities £H(t), &{t), and f(r) were defined in (|44| - (|46|) . 

We have now successfully controlled the required horizontal energy quantities. 
From Section EOl it remains only to control the diameters and the injectivity radii 
of the E T 's. First of all, Propositions [T7l and [521 imply uniform bounds for 

IWIi-(E T )i 70<T<tl. 

As a result, combining the above with Proposition[]j)l we obtain the desired uniform 
control for the injectivity radii of the E T 's. E3 

The diameters of the E T 's can be controlled by even cruder means using Propo- 
sition [131 For example, let p, q G S r , where To < r < t\. Consider the normal 
transports po , <7o of p, q to E To , and fix a curve ao in E To from po to go whose length 
is controlled by the diameter of E To . By using the transported coordinate systems 
of Proposition [13] and the uniform cllipticity property of (|64|) , then we can bound 
the length of the normal transport a of ao to E r . This bound depends only on the 
fundamental constants, including the diameter of E To . This controls the diameters 
of the E T 's and finally completes the proof of Theorem [6] 
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